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We perform observational tests of statistical isotropy using data from large-scale-structure surveys 
spanning a wide range of wavelengths. Using data from 2MASS, 2MRS, and NVSS galaxies, and 
BATSE gamma-ray bursts, we constrain the amplitude and direction of dipolar modulations in 
the number count of sources projected along the line of sight. We pay particular attention to 
the treatment of systematic errors and selection effects, and carefully distinguish between different 
sources of dipole signal previously considered in the literature. Dipole signals detected in these 
surveys are consistent with the standard, statistically isotropic expectation, except for the NVSS 
result, which is likely biased by remaining systematics in the data. We place constraints on the 
amplitude of any intrinsic dipole driven by novel physics in the early universe. 



I. INTRODUCTION 

The cosmological principle holds that the universe 
is homogeneous and isotropic on its largest observable 
scales. While the cosmological principle is a crucial in- 
gredient in obtaining many important results in quan- 
titative cosmology, there is no fundamental reason why 
our universe must obey it. The most general aim of the 
research we present here has been to directly test the cos- 
mological principle using data from recent astrophysical 
observations. 

Observations do bear out that the cosmological prin- 
ciple is a reasonable approximation for most purposes. 
However, few rigorous observational tests have been ap- 
plied to test homogeneity and isotropy. The work pre- 
sented here is directed toward performing tests of the sta- 
tistical isotropy of the universe using large-scale structure 
(LSS). The goal of this work is, fundamentally, to bring 
statistical isotropy more fully out of the realm of assump- 
tion and into the realm of observation. Any observed vio- 
lations of statistical isotropy could have far-reaching im- 
plications for our understanding of the universe's earliest 
moments, and violations of isotropy would also invalidate 
basic assumptions that serve as prerequisites to typical 
methods of data analysis in observational cosmology. 

To define statistical isotropy, consider a fluctuating 
field on the sky T(n) (the same arguments will apply 
for any field, including the cosmic microwave background 
(CMB) temperature field, the galaxy density field, etc.). 
The field is statistically isotropic if the two-point corre- 
lation function depends only on the separation between 
points: 

(T(n)T(n'))=C(n.n') (1) 

while in the case where statistical isotropy is violated, 
the right-hand side would read C(n, fi'). 

Alternatively, we could expand the field in spherical 
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harmonics T(n) = ^^T^ = a^^y£^(n), where 

i^O i^O m=-i 

the a£rn are the coefficients of the expansion. If statisti- 
cal isotropy is assumed, the angular power spectrum Ce 
may be defined via (a|^a£/^/) = CiSufSmrn' ^ where 
does not depend on the rotational degrees of freedom m. 
Hence, it is meaningful to calculate Q by averaging over 
the {2£ -\- 1) (on a full sky) "samples" corresponding to 
the {2i-\- 1) values of m for each i. However, without the 
assumption of statistical isotropy, we have in general 

i^im^^'m') = Clral'm'-) (2) 

which is much more difficult to measure than since we 
get only one sample of Cirni'm' for each (i^m^i^m^) in 
our universe. This is a big reason why statistical isotropy 
is such a crucial assumption in cosmology: it is much 
easier to work with Ci than the much more complicated 
quantity above. 



A. Point of Departure: Large-Scale Structure as 
Probe of CMB Anomalies 

While measurements of the angular power spectrum of 
the CMB by the Wilkinson Microwave Anisotropy Probe 
(WMAP) experiment have provided strong support for 
the inflationary Hot Big Bang model and allowed for 
unprecedentedly precise determination of cosmological 
parameters, these successes have come along with var- 
ious unexpected "anomalies" in the data (for reviews, 
see Copi et al. [1 and Bennett et al. [2 ). 

There are already tantalizing hints of violations of sta- 
tistical isotropy in CMB data; in fact, several observed 
anomalies, especially those indicating correlations be- 
tween patterns in the CMB of (apparently) cosmologi- 
cal origin (e.g., the quadrupole and octopole) and the 
geometry and motion of the Solar System, can be un- 
derstood as indications that statistical isotropy fails to 
hold. In addition, the southern ecliptic hemisphere has 
significantly more power than the northern ecliptic hemi- 
sphere on scales of about 3 degrees and larger (multipoles 
< 60) (Hoftuft et al. [3 ). This so-called "hemispherical 
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power anomaly" is a dipole modulation (see the begin- 
ning of Sec. [ni] for precise discussion of what this entails 
mathematically) of CMB power (direction in Galactic co- 
ordinates is (/,6) = (224°, —22°)) that is completely dis- 
tinct from the measured CMB dipole Ci (which is due to 
our motion with respect to the CMB rest frame), and is 
yet another possible indication of breaking of statistical 
isotropy in the CMB^. Finally, the low power at large 
angles in C{6) may itself be an indication that statistical 
isotropy is violated: it seems that a conspiracy of low-£ 
multipoles in the angular power spectrum is responsi- 
ble for creating the suppressed C {0) [5 , and correlations 
between different multipoles could be the result of a lack 
of statistical isotropy. 



B. Goals of This Work 

We take these considerations as general motivation for 
the work performed in the rest of this paper, which con- 
cerns testing statistical isotropy with large-scale struc- 
ture, and will analyze the specific issue of dipole patterns 
in various surveys being used as probes of the statistical 
isotropy of the universe. This work also complements the 
(surprisingly few!) studies of tests of statistical isotropy 
with current or future LSS observations [6HTQ]. 

The paper is organized as follows. In Sec. [llj we point 
out the various effects that are expected to contribute 
to a dipole signal in a large-scale-structure survey. In 
Sec. |III[ we outline the formalism used in this paper 
to detect such dipole signals. We then apply that for- 
malism to several surveys: the 2MASS Redshift Survey 
(Sec.|lV|), the Two Micron All-Sky Survey (2MASS) as a 
whole (Sec. [V|), the Burst And Transient Source Exper- 
iment (B ATSE ) of the Compton Gamma-Ray Observa- 
tory (Sec. [VIl| , and the NRAO VLA Sky Survey (NVSS; 
Sec. |VIIl[ ). In Sec. |VI[ we also examine issues surrounding 
searching for dipole signals in high-redshift objects and 
review work that has searched for dipoles in X-ray sur- 
veys. We conclude in Sec. |lXj summarizing our results 
in Table I Villi in that section. 



there are local anisotropics since the universe is not ho- 
mogeneous and isotropic except on its very largest scales, 
and (b) the Earth has a total motion relative to the large- 
scale-structure rest frame that is the sum of several vec- 
tor contributions (Earth moves around the Sun, the Sun 
moves around the center of the Milky Way, the Milky 
Way moves with respect to the Local Group barycen- 
ter, and the Local Group barycenter moves with respect 
to the structure around it and, ultimately, the large- 
scale-structure rest frame). That motion produces dipole 
anisotropy due to two effects, the Doppler effect and rela- 
tivistic aberration of angles (Itoh et al. [11]; see Sec. II C 
for mathematical details). The Doppler effect is relevant 
because it changes how magnitude varies with frequency, 
and since LSS surveys invariably operate within limited 
frequency ranges, the Doppler effect may shift certain 
objects into or out of a magnitude-limited sample. Since 
frequencies will increase in the direction of motion and 
decrease in the opposite direction, this produces a small 
dipole in the number of objects detected. Meanwhile, 
relativistic aberration causes the measured positions of 
galaxies to be displaced toward our direction of motion. 
This effect is on the order of v/c ^ 10~^, relevant for our 
purposes. 

We expect that as we go from smaller to larger sur- 
vey volumes, the measured value of the dipole amplitude 
should converge to that of the CMB dipole. This is be- 
cause with larger survey volumes, we average over more 
and more structure, and the universe approaches homo- 
geneity and isotropy. Any dipole left over should be due 
only to our motion, a kinematic dipole (with amplitude 
on the order of 10~^) just as in the CMB. There are sev- 
eral reasons why the dipole might not converge to that 
of the CMB: 

1. the rest frame of the CMB may not be the same 
as the rest frame of the LSS: novel horizon-scale 
physics (explored below) could induce a relative ve- 
locity between the CMB and LSS, so that galaxies 
would have a nonzero average streaming velocity 
with respect to the CMB rest frame; 



II. TYPES OF DIPOLES: SPECIFIC 
THEORETICAL CONSIDERATIONS 

It is completely expected that a dipole will be present 
in any survey of objects that trace large-scale structure. 
Both of the following effects contribute to the dipole: (a) 



2. there is also the possibility that there is genuinely 
more mass (and therefore more galaxies/objects 
that trace the mass distribution) in one direction, 
corresponding to modulation of primordial curva- 
ture perturbations due to the physics of infiation. 
For example, isocurvature perturbations can pro- 
duce such an effect (and explain the CMB hemi- 
spherical power anomaly) p!2UT5] . 



However, see Hanson and Lewis 4 for the updated result that 
the significance of the hemispherical power anomaly decreases 
when smaller scales are taken into account: previous analyses had 
been done for i < 60; when analysis is extended out to ^ ~ 500, 
the effect becomes less than 3cr-anomalous. The fact that greater 
resolution reduces the significance of the signal calls into question 
the authenticity of the signal as a genuine cosmo logical effect. 



Careful measurement of dipoles in various surveys, such 
as those we perform here, help zone in on these pos- 
sibilities, which correspond to a violation of statistical 
isotropy. 

The rest of this section will fiesh out the details of the 
discussion in the preceding two paragraphs. 



3 



A. Flux- Weighted Dipole vs. 2D-Projected Dipole 

One very commonly computed type of dipole is not, 
strictly speaking, a dipole at all, but is frequently referred 
to as such. This is the flux-weighted "dipole," where 
instead of calculating a (genuine) dipole based only on 
the two-dimensional projected positions of objects on the 
sky, some radial information is preserved by weighting 
each object by the flux we receive from it. We follow 
Bilicki et al. [16 in the following explanation of how the 
flux-weighted dipole is calculated. 

The flux-weighted dipole, as typically computed, is a 
measure of the acceleration due to gravity on the Local 
Group. In linear theory, the peculiar velocity v(r) at po- 
sition r is proportional to the peculiar acceleration vector 
g(r) induced by the matter distribution around position 
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Here Hq = 100 /ikm/s/Mpc is the Hubble constant, 
is the matter density divided by the critical density, and 
fiVtu) = {d\n D / d\n z)\z=o (where D is the growth fac- 
tor). In the ACDM model, the factor /(I^m) ^ 
and is only weakly dependent on the cosmological con- 
stant (Lahav et al. [18^). The acceleration vector itself is 
given by 



g(r) = Gp / (5m (r') 



3^/ 



(4) 



where Sm{'^) = [pm(i*) — p] /p is the density contrast of 
the mass perturbations at the point r; note that Pg is 
the mass density of galaxies, and b is the bias factor that 
relates mass density of galaxies to that of matter, b = 
Pg/ pM (assuming constant, scale- and time-independent 
bias). The bias is usually packaged with the factor /{^m) 
into the parameter (3 = f{^M)/b. Comparing Eqs. (|3| 
and (pi), we get the proportionality valid in linear theory: 



V (X /3g . 



(5) 



Comparison of the peculiar velocity v (determined from 
either the LSS surveys or using the CMB dipole) and 
acceleration g of the Local Group serves as a tool to esti- 
mate the P parameter. Independent knowledge of biasing 
allows one to estimate the cosmological matter density 
The programme of measuring the matter density 
in this way has been ongoing for over three decades (e.g. 

nMs,). 

The challenging part in this procedure is evaluating 
the acceleration in Eq. Q. This is where the flux- 
weighted-dipole approach comes in, where position vec- 
tors of the objects in a survey are weighted according 
to their fluxes (used as a rough proxy for mass, since 
both gravity and flux go as 1/r^), and added together. 
By preserving some radial information in this way, the 
flux- weighted dipole allows one to obtain a measure of 



the direction and strength of the acceleration of the Lo- 
cal Group due to the Newtonian gravitational attraction 
from objects in the survey. The flux- weighted dipole from 
2MASS has previously been found to be in the direction 
{l,b) = (264.5°, 43.5°) ± (2.0°, 4.0°) (Mailer et al. [29 ; 
note the rather serious discrepancy between the pub- 
lished result and the arXiv version, and see Bilicki et al. 
[16] for detailed discussion of convergence of the flux- 
weighted dipole and possible shortcomings of Mailer et 
al.); the flux- weighted dipole from 2MRS is in the direc- 
tion {l,b) = (251°, 38°) (Erdogdu et al. Hi]) in the Local 
Group frame, or (245°, 39°) in the CMB frame. 

Just to be clear, we do not consider further the flux- 
weighted dipole in this work; we henceforth study the 
2D-projected dipole described below. 



B. 2D-Projected Dipole: Local-Structure Dipole 

For the rest of this section, we focus on what we term 
the 2D-projected dipole, which is the quantity that is 
usually indicated by the isolated use of the word "dipole." 
This quantity relies on objects at any given redshift being 
projected on the celestial sphere (hence "2D-projected") 
with no weighting scheme. 

For a survey with very large (hundreds of Mpc- to Gpc- 
scale) volume, the universe is at least close to homoge- 
neous and isotropic on the scales relevant for the sur- 
vey. We naturally expect that any dipole signal in such 
a large-volume survey will be strongly suppressed. How- 
ever, on much smaller scales, where the universe is not 
at all homogeneous and isotropic, dipole signals should 
naturally emerge in any survey of objects that trace large- 
scale structure at all, and certainly in any galaxy survey. 
To take a particularly simple example, there is a large 
dipole in the galaxy distribution if we survey only ob- 
jects within the Local Group. 

But even if a survey encompasses structure on scales 
of tens of Mpc, we fully expect that given the non- 
uniformity of nearby structure, there should be a dipole 
component in the pattern of galaxies observed on the sky. 
This dipole component will be strongest for the smallest 
surveys, and should die away monotonically (at least in 
statistical average) for larger and larger surveys. The ef- 
fect turns out to be on the order of 10~^ for scales of 
tens to a couple hundred of Mpc, meaning that the fluc- 
tuations (contributing to the dipole) in the number of 
galaxies, as a function of position on the sky, are on the 
order of 1/10 the size of the mean number of galaxies 
across the sky. 

More rigorously, we make predictions for the full an- 
gular power spectrum of large-scale structure as a 
function of maximum redshift of a survey. The angular 
power spectrum of density fluctuations of halos is usu- 
ally expressed within the Limber approximation, where 
the contribution of modes parallel to the line of sight is 
ignored. In this approximation, the angular power spec- 
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trum is given by 



Assuming n(z) peaks at z=0.6 and bias=1.0 



Cf,= 
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where A^(/c) = /c^P(/c)/(27r^) is the dimensionless power 
spectrum, r{z) is the comoving distance, and H{z) is the 
Hubble parameter. The weight W{z) is given by 



W{z) 



b{z) 



dN 
dz 



dN 
dz 



(7) 



dz 



where z^in and z^ax 



are the lower and upper end of the 
redshift range, and dN/dz{z) is the number of galaxies 
per unit redshift.^ We adopt the tabulated values, or 
else functional form, of N{z) directly from the respective 
surveys that we study. 

However, we are interested in the dipole i = 1 where 
the Limber approximation is not accurate anymore (it is 
accurate at ^ > 10); see Fig. [l] Therefore, we adopt the 
exact expression for the power spectrum; using notation 
from (e.g.) Hearin et al. [32 , this is 

d = Att d\nkA'^{k,z = 0)l'^{k) (8) 
Jo 



I{k) = / dzW{z)D{z)j^{kx{z)) (9) 
Jo 

where x(^) is the radial distance, and x(^) = ^(^) 
a flat universe, the case that we consider. Here D{z) 
is the linear growth function of density fluctuations, so 
that 8{z) = D{z)5{0), where D{0) = 1. Note that, over 
the shallow range for 2MASS we can assume that b{z) is 
constant, and factor it outside of Eqs. Q and ([9|, but 
over the much deeper range for NVSS the bias may vary 
with redshift, and we adopt the expression for W{z) from 
Ho et al. [31 that implicitly integrates bias and number 
density as per Eq. This is explained in detail in 

Sec. lVniBl 

To produce the flducial theoretical predictions, we con- 
sider the standard cosmological model with the follow- 
ing parameter values: matter density relative to criti- 
cal Qm = 0.25, equation of state parameter w = —1^ 
spectral index n = 0.96, and amplitude of the matter 
power spectrum \nA where A = 2.3 x 10~^ (correspond- 
ing to cTg = 0.8) deflned at scale k = 0.002 Mpc"\ The 
power spectrum A^(A:,z) = k^ P(k) / {2ti'^) is calculated 
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FIG. 1: A plot of the angular power spectrum Cg predicted 
for a galaxy survey with a peak in the galaxy redshift distri- 
bution at z — 0.6 and the given maximum depth Zmax- In 
the remainder of this paper we will focus attention on the 
dipole, i — \. For this particular redshift distribution, the 
local-structure dipole becomes subdominant to the kinematic 
dipole for around Zmax 1.0. For higher Zmax, we should get 
convergence to the kinematic dipole plus any intrinsic dipole 
that might be present. The dotted curves correspond to the 
power spectrum within the Limber approximation (low ^) and 
assuming linearity (high ^); the solid curves correspond to the 
more accurate set of assumptions where the Limber approxi- 
mation is relaxed and nonlinearity is taken into account. 



using the transfer function output by CAMB. We do not 
vary the values of cosmological parameters, since they 
are measured to suflicient accuracy that any shifts in pre- 
dicted dipole amplitude that could occur due to realistic 
changes in cosmological parameters are tiny in compar- 
ison with cosmic variance given the flnite sky coverage 
and relative shallowness of the surveys we employ (as we 
have explicitly verifled). 

The power spectra are shown in Fig. [l] Note the sub- 
stantial, order-unity, difference between the exact and 
approximate (Limber) expressions at ^ = 1; in the re- 
mainder of this work we use the exact, double-integral 
expression. 



C. 2D-Projected Dipole: Kinematic Dipole 



Note that a sometimes- used alternative definition of n{z) refers 
to the spatial density of galaxies (e.g. Hu and Jain [30|); it is 
related to the quantity we use via dN j dz = n(z) Q r'^ (z) / H (z) , 
where Q is the solid angle spanned by the survey, and r and H 
are the comoving distance and Hubble parameter, respectively. 
Note also that our W{z) is equivalent to the quantity f(z) from 
Ho et al. ^31,. 



A dipole pattern may also arise due to motion of the 
Earth with respect to the astrophysical objects or struc- 
ture being measured. This is what produces the dipole in 
the cosmic microwave background, and it also contributes 
to the total dipole in a measurement of the large-scale 
structure. 
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TABLE L Motions that give rise to a kinematic dipole in the CMB and large-scale structure. 



Motion 


Approximate Speed (km/s) 


Direction 


Earth around Sun 


- 30 km/s 


annually varying 


Sun wrt Local Group 


- 306 km/s 


(/,6) = (99,-4)±(5,4) 


Local Group wrt CMB 


- 622 km/s 


{l,b) = (272,28) 


Overall CMB kinematic dipole 


- 370 km/s 


(/,6) = (264.4,48.4) zb (0.3,0.5) 



1. Kinematic Dipole in the CMB 

Probably the best-known dipole in all of cosmology is 
the dipole measured in the CMB temperature distribu- 
tion. This dipole, which has an amplitude on the order of 
10~^ times the amplitude of the CMB monopole, arises 
due to the motion of the Solar System with respect to 
the CMB rest frame. This motion is the vector sum of 
several different motions, summarized in Table |lj 

Values of the kinematic dipole in the CMB are cited 
with the contribution from the Earth's motion around 
the Sun subtracted out, so that the dipole is due only 
to the Sun's velocity with respect to the CMB (Kogut 
et al. [33^). The value of the Local Group's peculiar 
velocity with respect to the CMB is from Mailer et al. 
[29^ and was computed using the value of the Sun's ve- 
locity with respect to the Local Group in Courteau and 
van den Bergh [34]. (When the velocity of the Local 
Group with respect to the CMB rest frame is inferred 
from the measurement of the CMB dipole, the direction 
becomes (/,6) = (276,30) ± (3,3) [33]. Compare also 
values determined in Tully et al. [35] and Bilicki et al. 
[16].) The peculiar velocity predicted from linear-theory 
ACDM is ^ 470 km/s [36 . Note that the speed of the 
Sun with respect to the CMB rest frame would be consid- 
erably greater if not for the fact that the Sun's velocity 
vector with respect to the Local Group points in a direc- 
tion nearly opposite that of the Local Group's velocity 
vector with respect to the CMB. Also note that the domi- 
nant contribution to the Sun's motion with respect to the 
Local Group is the Sun's motion around the center of the 
Galaxy, which has speed ~ 220 km/s, and is composed of 
the Sun's motion with respect to the Local Standard of 
Rest and the LSR's motion with respect to the Galactic 
Center [111|34]. 

2. Kinematic Dipole in LSS 

The kinematic dipole in the CMB, which is due to 
the Sun's motion with respect to the CMB rest frame, 
is observed as a Doppler shifting of the CMB photons. 
The effect that gives rise to a kinematic dipole in the 
large-scale structure is not quite as direct. Rather, it 
includes contributions both from the Doppler effect and 
relativistic aberration. We derive the relevant equations 
in Appendix [B] 



D. 2D-Projected Dipole: Intrinsic Dipole 

In the CMB, the intrinsic dipole corresponding to adi- 
abatic perturbations is zero (Erickcek et al. [12 ). When 
we switch over from talking about the CMB to talking 
about large-scale structure, we expect that there may be 
an intrinsic dipole in the LSS. Below, we explore possible 
reasons why an intrinsic dipole might compete with or 
even (conceivably) dominate the LSS kinematic dipole. 

Erickcek et al. [12 propose a scenario in which the 
curvaton (particle mediating a scalar field that may gen- 
erate fluctuations during inflation without actually driv- 
ing inflation) has a large-scale spatial gradient, which in 
turn causes variation in the amplitude of the primordial 
curvature perturbations, modulating A^^ across the sky. 
Hirata [7 shows how this modulation due to isocurva- 
ture perturbations would transfer to the CMB and large- 
scale structure, in the latter case causing a dipolar vari- 
ation in the abundance of massive haloes (and objects 
that occupy them). This inflationary scenario is one sce- 
nario that invokes the physics of the early universe to 
explain why there might be an intrinsic dipole in the 
large-scale structure above and beyond what we naturally 
expect to be present from typical scale- invariant fluc- 
tuations/adiabatic perturbations laid down in the sim- 
plest inflationary scenarios. While the simplest curvaton- 
gradient model has been ruled out by Hirata's analysis 
of constraints on the dipole in SDSS quasars, and cor- 
responding constraints on dipolar modulation of the pri- 
mordial power spectrum,^ similar but more complicated 
scenarios are still possible. 

Note that Hirata's constraints on the primordial dipole 
amplitude using SDSS quasars are on the order of 2 x 
10~^, which corresponds to constraints on the amplitude 
of the dipole in the quasars themselves roughly an order 
of magnitude higher. Hence current constraints on this 
particular intrinsic-dipole scenario are not down to the 
level associated with the kinematic dipole, though this 
was not a problem in Hirata's analysis since that anal- 



^ According to Hirata, any smooth gradient in the amplitude of the 
primordial curvature perturbations is no more than 2.7 percent 
per present-day horizon radius (99 percent confidence); cf. the 
11 percent variation required in the Erickcek et al. model needed 
for consistency with the CMB hemispherical power anomaly. 
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ysis looked specifically for a dipole effect that would ac- 
company the primordial conditions needed to explain the 
CMB hemispherical power anomaly given the curvaton- 
gradient model, and that would have required a 10~^ 
dipole. 

Another possibility for generating an intrinsic dipole is 
that the CMB rest frame is not the same as the large- 
scale-structure rest frame. This happens, for example, in 
models with a large Gpc radius) underdense void, in 
which we are located close to the center. The observed 
bulk flow is then equal to the difference between the Hub- 
ble parameters inside and outside the void multiplied by 
our distance away from the center^ |37j. Requiring that 
the intrinsic dipole thus measured is consistent with ob- 
servations of the CMB dipole (v/c ^ O(10~^)) requires 
that we live very close (< 15Mpc) to the center of the 
void [38], making these models rather fine-tuned. 

An alternate mechanism for how a CMB-LSS rest- 
frame disagreement may arise is provided by the 
Grishchuk-Zel'dovich effect (Grishchuk and Zeldovich 
[39 ; Erickcek et al. [40 j see also Gunn [41 ). As Turner 
[42j explains, if inflation lasts only a little longer than 
necessary to solve the flatness and horizon problems, 
scales that were superhorizon-sized at the onset of in- 
flation - these scales cannot be affected by events dur- 
ing or later than the inflationary epoch, and thus con- 
tain imprints of the pre- inflationary universe - may not 
be much larger than our present horizon, and thus may 
have some effect in the current universe. In particular, 
he proposes that large density fluctuations with wave- 
lengths slightly larger than the Hubble radius (modes 
that are "just barely" superhorizon-sized) may exist, and 
would appear to us as a density gradient in a partic- 
ular direction. Such a density gradient could produce 
a "tilted universe": a universe in which all the matter 
within the Hubble volume gains a peculiar velocity due 
to the greater gravitational attraction from one "side" of 
the universe than the opposite. The effect is equivalent 
to saying that the rest frame of the CMB is not the same 
as the rest frame of large-scale structure: from the rest 
frame of the CMB, all matter would have a nonzero av- 
erage streaming velocity. This would produce an intrin- 
sic dipole due to relativistic aberration and the Doppler 
effect (or, equivalent ly, it would produce an additional 
kinematic dipole on top of that expected from analysis 
of the CMB) [IJ. In order for this "tilting" effect to be 
observable, isocurvature modes must be present, even in 
the presence of late-time acceleration [43 . 

Given the presence of isocurvature modes, the 
Grishchuk-Zel'dovich effect would also produce an (ad- 
ditional) intrinsic dipole due to the simple fact of the 



Therefore, the CMB frame and the LSS frame inside the void 
(but not the global LSS frame) are different in this case. Most 
measurements of bulk flows, being much shallower than the size 
of the void, would interpret this as a legitimate difference be- 
tween the CMB and LSS frames. 



superhorizon-scale density perturbation. As another ex- 
ample of physical mechanisms that would produce an in- 
trinsic dipole, Gordon et al. [44] examine a scenario in 
which there are spatial perturbations in the density of 
dark energy from a quintessence field: that is, a late-time 
effect produces horizon-scale fluctuations. More gener- 
ally, Gordon et al. examine a class of models in which 
the full fundamental theory is homogeneous and statis- 
tically isotropic, but statistical isotropy is broken from 
a given observer's position because of superhorizon-scale 
perturbations that appear as a gradient in density across 
the sky on the largest observable scales. Any theory that 
generates such a variation in density would give rise to 
what we have termed an intrinsic large-scale-structure 
dipole, and the appearance of the breaking of statistical 
isotropy. These density variations could, at least theo- 
retically, exist on essentially any order of magnitude in 
6p/p. 

There is some reason to take the idea of a tilted uni- 
verse seriously. Kashlinsky et al. [45 investigate the bulk 
motion of galaxies in the universe out ~ 300 Mpc/h and 
find, somewhat controversially (see, e.g., [46-48 ), that 
there is a coherent bulk flow in their sample. The evi- 
dence they develop for this claim comes from attempting 
to detect the kinetic Sunyaev-Zel'dovich effect by com- 
puting the dipole of the CMB temperature field eval- 
uated at the positions of galaxy clusters. This dipole, 
evaluated as it is in a small number of pixels, does not 
receive appreciable contributions from our own motion 
(i.e., from the CMB kinematic dipole due to the Sun's 
motion with respect to the CMB rest frame), but does 
receive contributions from instrument noise, the thermal 
SZ effect, the intrinsic CMB dipole, and foreground com- 
ponents. However, contributions other than the kine- 
matic SZ effect are, they argue, accounted for in their 
analysis, with the thermal SZ effect in particular can- 
celing out /integrating down when averaged over a large 
number of clusters. Their conclusion is that the dipole in 
CMB temperature evaluated at cluster positions is due 
to the kinetic SZ effect due to the bulk flow of the clus- 
ter sample. If this effect is authentic, then it fits well 
with the tilted-universe scenario: the bulk motion is de- 
tectable in large-scale structure but does not generate a 
primordial dipole CMB component. 

Other studies that use measurements of peculiar ve- 
locities of local neighborhood galaxies — typically de- 
termined by combining the measurements of their dis- 
tances and redshifts — find larger-than expected bulk 
flows [49','50 . These flows are of order 400 km/s and seem 
to be showing no signs of convergence out to the probed 
distance R ^ 60 /iMpc~^; they are estimated to be ^ 1% 
likely in the standard ACDM cosmological model for the 
given observed volume. The relation of these larger-than- 
expected bulk-flow measurements to findings by Kash- 
linsky et al. [45 is unclear at this time [50 , especially 
given that bulk flows inferred from distances obtained 
from type la supernovae indicate somewhat lower bulk 
flows that are therefore in better agreement with ACDM 
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[5TH53]. 

E. Types of Dipoles: Review 

In summary, when we observe some tracer of large- 
scale structure (galaxies, quasars, gamma-ray bursts, 
etc.), we may observe a dipole in counts. If the dipole we 
are observing is what we have called the 2D-projected 
dipole - that is, a dipole in surface density of the ob- 
ject - then contributions to that dipole may come from 
(1) the local-structure dipole, (2) the kinematic dipole 
(which is due to the Doppler effect and relativistic aber- 
ration), and (3) an intrinsic LSS dipole (which is really 
the z ^ 1 limit of the local-structure dipole). There are 
a couple observations to make about these effects: 

• In the limit of very high redshift, the local-structure 
dipole goes to amplitudes on the order of 10~^ (as 
we have explicitly verified), the kinematic dipole 
is expected to go to amplitude vjc 10~^ and 
align with the direction of the CMB dipole, and 
the intrinsic dipole could take on a wide variety of 
values depending on certain theoretical considera- 
tions. 

• For structures/galaxies at relatively small redshifts, 
the local-structure dipole amplitude >> the kine- 
matic dipole amplitude. However, even though the 
kinematic dipole is swamped by the local-structure 
dipole, we expect that these two dipoles should 
point in somewhere close to the same direction. 
While no particular level of agreement is guaran- 
teed, the fact remains that local structure is what 
accelerates us in the direction that the kinematic 
dipole points. This is why, in linear theory, the 
velocity of the Local Group is proportional to its 
acceleration due to gravity. However, since the 2D- 
projected dipole takes no radial information into 
account, it is not a true measure of gravitational 
attraction or acceleration, but only a partially reli- 
able proxy. 

In Table [HI we show the types of dipoles considered in 
this work, and the section in which they are introduced. 

III. FORMALISM FOR DETECTING DIPOLES 

Some forays have already been made into tests of sta- 
tistical isotropy, and dipoles in particular, using measure- 
ments of large-scale structure. Many estimators for the 
dipole have been employed, some of which do better jobs 
than others at naturally incorporating sky cuts, allowing 
for systematic effects to be accounted for, etc. Here we 
adopt the estimator used by Hirata [7 to test the WMAP 
hemispherical power anomaly using quasars detected by 
the Sloan Digital Sky Survey. This is the best unbiased 
estimator for determining the amplitude and direction of 



a dipole in counts of objects on the sky under conditions 
of cut skies and in the presence of systematics. We now 
describe this formalism. 



A. Obtaining the Direction of the Dipole 

Consider a dipolar modulation on the sky with some 
amplitude A in a (unit) direction d. We may write the 
observed density field N of the objects in question as a 
function of direction n as 

N(xi) = [1 + A(d • n)]7V + e(n) (10) 

where N is the intrinsic statistically isotropic field and e 
combines random and instrumental noise. If we momen- 
tarily drop the e term, we can write SN/N = A{d • n). 
Reinstating a term corresponding to systematic errors, 
the fluctuations in density as a function of direction can 
be written as the sum of contributions from a dipole, 
fluctuations due to systematics, and a mean offset [7] 

^(fi) = Ad • n + ^ kMn) + C. (11) 

i 

Here ti{n) are possible systematics templates in the sky 
map (such as an extinction map), the coefficients ki give 
the amplitudes of the contributions of these systemat- 
ics to the observed density field, and the presence of the 
monopole term, C, allows us to account for covariance 
between the monopole and other estimated parameters, 
especially covariance between the monopole and any sys- 
tematic templates. 

It is then straightforward to write down the best linear 
unbiased estimator of the combination (d, /c^, C) with 
corresponding errors. The procedure is as follows: First, 
we rewrite the above equation as 

SN 

^(n)=x-T(n) (12) 

where x = {dx,dy,dz,ki, ...,kN,C), T(n) = 
{nx,ny,nz,ti{n), ...,tN{fi), 1), and + + = 1. 
The best linear unbiased estimator of x is 

^ = F-^9 (13) 

where the components of the vector g are 

gi = j Ti{n)dN''{h)d^n (14) 

and the Fisher matrix F is given by 

F^J =N'' J T,{h)Tj{h)d^h, (15) 

where = dN/dQ is the number of galaxies per stera- 
dian. To actually compute these quantities with dis- 
cretized data, it is convenient to work with a data map 
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TABLE II: Table reviewing the sources of dipole signal in a large-scale- structure survey. 



Local-Structure Dipole 



Kinematic Dipole 



Intrinsic Dipole 



Typical size 



0.1-10 depending on Zn 



10" 



>io- 



Flux- Weighted 
method probes: 



Local Group's^ 
acceleration: Sec. 



II A 



N/A 



N/A 



2D-Projected 
method probes: 



small-scale departures from 
statistical isotropy; Sec. II B 



our motion; 



Sec. lie 



many possible the oretic al 
explanations; Sec. II D 



and a random map, the latter of which is simply a set of 
randomly chosen directions/points hji on the unit sphere: 



9i 



D 



Nr 



R 



Nd 

Nr 



(16) 



(17) 



where Njj and Nji represent galaxy counts rather than 
the number of galaxies per steradian as in the continuous 
case. 

Note that the component of g corresponding to the 
monopole term in Equation ( 11 ), which we will refer to as 



gcj must be zero, even if the sky is cut. This can be seen 
in the analytic formula for g by noting that we are inte- 
grating fluctuations relative to the mean, where the mean 
is determined from whatever portion of the sky is being 
integrated over. In the formulation where we discretize 
the celestial sphere, gc = Nd - {Nd/Nr){Nr) = 0; 
gc represents the monopole of the fluctuations from the 
mean on the cut sky, which must be zero. Hence, the 



only way the monopole term C in Equation (11) can be 



nonzero is by picking up on the covariances between vari- 
ables. 

To show explicitly how we calculate the Fisher matrix 
Fij in the discrete formalism, we take F^z as an example: 



Nd 
Nn 



i=i y^^R/ 



(18) 



where we have used z to designate the z-coordinate of the 
vector pointing to the center of the pixel in which count 
i is found. Since (2;^) = 1/3 over the entire sphere, we 
have Fzz = N^/^ (for the entire celestial sphere) in the 
limit of sufficiently large number of counts in the random 
map to have suppressed Poisson noise. 



variance matrix is CoY{pi,pj) = {F~^)ij. By the Cramer- 
Rao inequality, the best-case marginalized errors on the 
parameters are 



CTmargfe) = ^/{F^^i; 



(19) 



inverting F automatically mixes all the elements together 
and takes into account how they covary. Meanwhile, the 
best-case unmarginalized errors are 



r{pi) = l/\/Fii. 



(20) 



Note that the errors on our estimates of the dipole are 
based on the shape of the sky cut, the input systematic 
templates, and the number of data points N]j. 

As a side note, the correlation between parameters pi 
and pj is 



Fr 



(F-^FrM 

V II 11 ) 



(21) 



C. Estimating the Amplitude of the Dipole 

In this formalism, we need only to acquire the com- 
ponents of the dipole (dx^dy^dz\ and the associated er- 
rors [px^^yy^^yz)' Combining the components by squar- 
ing, summing, and taking the square root of the sum 
would create a biased estimator of the dipole amplitude 
A, so we never do this. Instead, once we have the best- 
fit dipole dbest = ^d, we can construct a marginalized 
likelihood function for the amplitude A [7 : 



-(An- dbest )Cov ^(Afi- dbest) 



(22) 

where Sxv indicates integration over all possible direc- 
tions on the sphere.^ In this equation, for a given am- 
plitude A, we take the best-fit dipole dbest as a given. 



B. Errors in Estimated Dipole Direction 



The matrix Fij is the Fisher matrix for the full param- 
eter set Pi = {dx^ dy^ dz^ /ci, /cat, C}, and hence the co- 



^ Even though c?best itself is a biased estimator of the amplitude 
A, our likelihood as written in Eq. (|22| returns an unbiased es- 
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and then compare each direction n (with the given am- 
phtude) with the best-fit dipole. The likelihood function 
is a Gaussian in A for fixed direction n, by construction, 
but may not be Gaussian when marginalized over direc- 
tion. That marginalization occurs in the equation above 
when we integrate over all fi, giving us the likelihood of a 
particular amplitude A marginalized over all directions, 
given the best-fit dipole dbest "selected" by the data. 
Posterior analysis will then show where 95 percent of the 
weight lies. 

Given that we ultimately work discretely, with a ce- 
lestial sphere that is pixellized using HEALPix the 
likelihood turns into a sum over pixels: 



For the purposes of order-of- magnitude calculations, the 
rule of the thumb is Ci A^. 



C{A) (X ^ exp 



-(An - dbest) Gov ^(An- dbest) 



E. Commentary on the Formalism 

It is straightforward to show that this estimator is ei- 
ther precisely or approximately equivalent to similar esti- 
mators used by recent authors. Nevertheless the formal- 
ism of Hirata that we presently use has several practi- 
cal advantages. First, the real-space estimator employed 
here is more convenient to implement than multipole- 
space estimators employed in previous analyses (e.g., 
AArea'^^^ et al. [55 , Blake and Wall [56 , Baleisis et al. [57] , 
etc.). Many analyses use pseudo-C^ to deal with sky cuts, 



(23) 

The factor AArea will come out of the summation since 
all pixel areas are equal in HEALPix and all that matters 
is the ratio of likelihoods rather than the absolute values 
of likelihoods, so we literally sum over all the pixels in 
order to get the marginalized C{A). We drop the pref- 
actor on the likelihood that includes covariance since the 
covariance does not depend on parameters. 

Finding the likelihood distribution as a function of di- 
rection, >C(n), follows from an exactly analogous proce- 
dure, but we sum over all possible amplitudes associated 
with a given pixel rather than over all possible pixels 
associated with a given amplitude. 



while sky cuts are very straightforward to deal with in 
this formalism (see Sec. HIE for further details). Finally, 
estimating the coefficients ki allows one to very natu- 
rally incorporate any systematics templates one suspects 
might be relevant and ensure that they do not interfere 
with estimation of the dipole. This form of component 
separation allows one to isolate the different contribu- 
tions to the observed ffuctuations in counts, and separate 
those contributions into actual dipole plus systematic ef- 
fects. Any pattern put into this formalism as a systematic 
template will be marginalized over in the determination 
of dipole amplitude A and direction d. 

Moreover, this formalism has several other noteworthy 
features: 



D. Converting From Dipole Amplitude A to 
Angular Power Spectrum Ci 

Here we show that there is a simple relationship be- 
tween the dipole power Ci, the £ = 1 mode of the an- 
gular power spectrum familiar from several areas of cos- 
mology, and the amplitude A of the dipole computed 
above. Without loss of generality, we assume that the 
dipole points in the positive z direction and write the 
ffuctuation in counts in two ways: 



SN 



(h) = Ah- d = A cose = aioYio(n) (24) 



where Yio = \/3/(47r) cos^. Therefore, aio = \/4:7t/3A. 
Using the usual relationship Ci = ^ \a£m\'^ /{2£ + 1), 



the power spectrum Ci contribution is then given by 



(25) 



timate. The reason is that the hkehhood effectively compares 
theory to measurement in each component (x, y, or z) separately, 
as can be seen by diagonalizing and rewriting the exponent in 
that expression, and the components themselves are unbiased. 
Numerical evidence that the estimate of A is entirely unbiased 
is shown in Fig. I20| in Appendix Ib| 



• It allows very naturally for arbitrary sky cuts: all 
that is necessary is to remove pixels from both the 
data map and the random map when performing 
the dipole analysis. When the sky is cut, the dipole 
becomes coupled to other multipoles, and the errors 
on the detection of dx, dy, and dz derived from the 
Fisher matrix correspondingly increase to account 
for this. 



• It allows for straightforward incorporation of arbi- 
trary pixellization. The scale of the pixellization 
should not matter, because Poisson noise is on the 
scale of the pixellization, which is much smaller 
than the scale of the dipole, and Poisson noise in 
larger pixels means a smaller effect (goes as 1/ Vn)^ 
so the Poisson noise cancels out. However, differ- 
ent pixellization schemes can affect the size, shape, 
and nature of a mask if there is any sky cut, and 
when discrepancies appear between dipole results 
using different pixellization schemes, it is virtually 
always traceable to this. 

• It allows for the possibility of A > 1. This 
may seem counterintuitive since it implies negative 
counts in some pixels. However, even though it is 
true that counts cannot go negative in real data, 
it is still possible for a model in which some pixels 
have negative counts to be the best fit to the data. 
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Note that the level of degeneracy between the system- 
atic template ti{n) and the dipole d depends on struc- 
ture of the former, and on the relative orientation of the 
two. We have performed extensive tests to verify and 
build our intuition about this degeneracy. For example, 
if a pure-dipole map is created with the dipole pointing 
in the z-direction, and this map is used as a template 
t(n), the effect is that the detection of the dipole in the 
z-direction, d^, gives an unreliable number, while the er- 
ror bar in the z-direction, a^, blows up to a much larger 
number than d^, so that {S/N)^ = dz/cTz ^ 1- Also, the 
correlation pzk between dz and the template coefficient 
k becomes 1.0. In other words, a dipole template in a 
given direction takes out any component of the dipole 
detection in that direction. 

Finally, we performed tests with varying sky coverage 
and survey depth, and the number of (mock) galaxies 
available to verify that the input dipole is successfully 
recovered within the estimator's reported error. 

We proceed to apply this estimator to data from sev- 
eral surveys. We select surveys with very wide sky cover- 
age, ideally almost full-sky coverage, because for a fixed 
depth, the number of modes available scales as the frac- 
tion of the sky covered, /sky This is useful especially for 
beating down cosmic variance in theoretical predictions, 
which as we will see, is the dominant source of uncer- 
tainty in our comparisons of observations with theory 
(observational results tend to be much more tightly con- 
strained than theoretical predictions, since we are work- 
ing at very low £). 

We also find that higher multipoles do not contribute 
appreciably to the recovered signal, or even strongly af- 
fect the error bars on the dipole signal, as long as more 
than roughly half the sky is probed in the survey. See 
Appendix [B] for details of how we use maps of ^ = 2 
and i = 3 modes as systematics templates to detect the 
presence of coupling among multipoles. 

Previous research on dipoles in similar surveys will be 
profiled as different types of surveys are brought up. 



IV. DIPOLE IN 2MRS 
A. Introduction to Dipole Signals in 2MASS 

The Two Micron All Sky Survey (2MASS), which im- 
aged 99.998 percent of the celestial sphere (Skrutskie 
et al. [58]), provides an excellent starting point in testing 
for dipoles in tracers of large-scale structure. This sur- 
vey includes two main catalogs, the point-source catalog 
(PSC) and extended-source catalog (XSC). The latter 
is of interest here, since it includes roughly 1.6 million 
sources, nearly all of which are extragalactic. 

2MASS used two 1.3-m equatorial Cassegrain tele- 
scopes, one in the Northern Hemisphere and one in 
the Southern Hemisphere (Mt. Hopkins, Arizona; Cerro 
Tololo, Chile), to observe in the J, and Kg bands, 
corresponding to wavelengths of 1.25, 1.65, and 2.16 /im. 



respectively. The XSC contains sources that are ex- 
tended with respect to the instantaneous point-spread 
function [58], including galaxies and Galactic nebulae. 
The S/N=10 sensitivity limits are met by sources as 
bright or brighter than 15.0, 14.3, and 13.5 mag in the J, 
and Ks bands, respectively, and (very importantly for 
our dipole-related considerations) exhibit a mean color 
difference of less than 0.01 mag between hemispheres, 
meaning that the photometry is highly uniform between 
hemispheres. The reliability (corresponding to the ratio 
of the number of genuine extended sources to the to- 
tal number of sources, spurious or genuinely extended, 
in the dataset) of the XSC is greater than 99 percent 
for Galactic latitude \b\ > 20°. Some extended sources 
in the catalog are not extragalactic, though these can be 
easily removed with the right color cuts (as detailed later 
in Sec.[V|). 

A small subset of the 1.6 million extended sources in 
the 2MASS XSC were assigned redshifts in the 2MASS 
Redshift Survey (2MRS), a catalog which includes po- 
sition and redshift information for over 40,000 galaxies 
present in the original 2MASS sample. In this section 
and the next, we apply the dipole-detecting formalism 



outlined in Sec. [Hi] to the entire 2MRS catalog, as well 
as to appropriate subsets of the 2MASS XSC sources. 

Erdogdu et al. [17 , Mailer et al. [29 , and Bilicki 
et al. [T6^, for example, have calculated the flux- weighted 
dipoles (see Sec. |IIA| ) for 2MRS (a 23,000-galaxy subset 
thereof, actually, with Ks < 11.25 - a preliminary version 
of the catalog) and 2MASS, respectively. This stands 
in a longer tradition of attempting to calculate flux- 
weighted dipoles from near-infrared surveys, since near- 
infrared light closely traces the mass distribution of large- 
scale structure. For instance, Rowan-Robinson et al. [59] 
calculate a flux- weighted dipole from the IRAS PSCz 
Redshift Survey, which had redshifts for over 15,000 
IRAS galaxies (at 60 /im; cf. the wavelengths of 2MASS, 
over an order of magnitude shorter). The IRAS PSCz 
(^max ^ 0.1), 2MRS (^max ^ 0.1), and 2MASS XSC 
{z > 0.07) studies all find tolerably small discrepancies 
between the direction of the flux-weighted dipole (and 
thus the acceleration of the Local Group) and the CMB 
velocity dipole that partially results from that accelera- 
tion (velocity of the Local Group being proportional to 
acceleration of the Local Group in linear theory). (As 
noted before, the motion of the Sun with respect to the 
Local Group also contributes to the kinematic dipole, but 
the direction is nearly opposite the direction of the Local 
Group's motion with respect to the CMB rest frame, and 
hence changes the magnitude of the velocity vector but 
does not substantially change its direction.) 

The number-weighted dipole of Erdogdu et al. [17] 
comes closer than the flux-weighted dipole to mimick- 
ing the quantity that we calculate here, but the number- 
weighted dipole, like the closely related flux- weighted 
dipole, is another quantity that seeks to measure the ac- 
celeration of the Local Group due to surrounding large- 
scale structure, but instead of using flux as a proxy for 
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mass as in the flux- weighted dipole, the number density 
of galaxies in a given direction on the sky is assumed to 
serve as a good proxy for mass. Our goals and aims, as 
well as the precise quantity we calculate, are different: 
We seek to measure not an acceleration, but rather the 
simple 2D-projected dipole, which in the case of the rel- 
atively nearby survey 2MASS (and 2MRS, as a subset of 
2MASS), is dominated completely by t he co ntributions 
from the local-structure dipole (see Sec. JIB). 



2MRS Galaxies 



The entire power spectrum of 2MASS has been calcu- 
lated by Frith et al. [55 , which means that at least one 
measure of the 2D-projected dipole that we explore here 
has already been obtained. We compare our results to 
this previous result later in this section. However, we do 
not regard our result as a simple replication of the pre- 
vious result (and, in fact, we note substantial disagree- 
ment): we compute not just the amplitude of the dipole, 
but also its direction; we account for systematics in a 
direct and natural way; and we place the 2MASS dipole 
into a larger context of exploring the various contribu- 
tions to dipoles, and testing observational results against 
theoretical predictions, in a wide variety of surveys. 




0.04 0.06 

Redshift z 



B. 2MRS Profiled 

We begin with the 2MASS Redshift Survey, the dens- 
est all-sky redshift survey to date. The 2MRS team 
(Huchra et al. [60]) measured redshifts of 43,533 bright 
\Ks < 11.75) sources with ^(5-1/) < 1 mag and \b\ > 5° 
(for 30° < / < 330°; |6| > 8° otherwise). Sources were 
carefully screened to ensure that all were genuinely extra- 
galactic sources and do not have compromised photome- 
try. As explained below, we err on the conservative side 
and make a symmetric cut at \b\ < 8° for all Galactic lon- 
gitudes which eliminates roughly 1700 of the galaxies 
in the survey. Nearly all 2MRS galaxies are within the 
range < z < 0.1. Previous tests have worked with the 
flux- weighted dipole in 2MRS [17 , but have not explored 
the various contributions to the 2D-projected dipole. 

The results for the dipole amplitude in this survey are 
strongly expected to agree with theoretical predictions. 
Given the relatively small volumes surveyed (and the fact 
that we are dealing with the very large-scale £ = 1 mode), 
cosmic variance is very large, and so a statistically signif- 
icant discrepancy between theory and observation would 
require highly anomalous disagreement between the two. 
We flnd, even with relatively cursory checks, that there 
are no serious discrepancies between theory and observa- 
tion for 2MRS. Nevertheless, in the coming sections, we 
proflle the various tests performed on the data, and the 
results for both dipole direction and amplitude, in the in- 
terests of presenting results for this survey as something 
of a model (in addition to being an important test in its 
own right): this is a dataset with well-controlled system- 
atics and very little chance of giving anomalous results, 
where we can test out several different types of system- 
atic checks, to gain intuition for what results should look 



FIG. 2: Top panel: All sources in the 2MASS Redshift Survey 
that escape the \b\ < 8° cut. The mean redshift in the survey 
is approximately z = 0.028. Even by eye, it is clear that the 
dipole due to local structure has not died away at these scales. 
In particular, the supergalactic plane is still fairly clearly vis- 
ible in the data (see, e.g., Mailer et al. |29]). (Note that 
the dynamic range of this plot has been limited so that struc- 
tures outside the supergalactic plane are also visible.) Bottom 
panel: A plot of the radial distribution of 2MRS galaxies. The 
data are put in redshift bins [0, 0.01), [0.01, 0.02), [0.09, 
0.10), where the plot shows the number of galaxies in each 
bin as a function of bin center. 



like when we perform similar tests on surveys at higher 
redshifts and/or with less well-understood systematics. 

This survey is essentially ready to be analyzed 
"straight out of the box," meaning that major systematic 
errors have already been addressed (especially Galactic 
extinction), and we already have a sample of extragalac- 
tic sources with uniform sky coverage outside the Galac- 
tic plane. (The latter is important since a lack of uni- 
form completeness across the sky could, if not properly 
accounted for, mimic the effect of a dipole.) More careful 
attention must be paid to these matters in the 2MASS 
sample as a whole, and in other surveys, but 2MRS re- 
quires only that we cut out pixels (in both the data map 
and the random map to which it is compared; see Sec. Ill) 
within 8 degrees of the Galactic equator, \b\ < 8°. Note 



that pixels whose centers are above 8 degrees, and thus 
escape a straightforward cut of pixels with centers below 
8 degrees, still may have area below 8 degrees, especially 
if the pixelization is coarse, as in the cases of HEALPix 
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Although 2MRS should be relatively systematic- free, 
we proceed to perform straightforward tests for two types 
of systematic effects: Galactic extinction as characterized 
by the maps of Schlegel et al. [61 , and star-galaxy con- 
fusion or other systematics that vary as a function of 
Galactic latitude b. 



D. Comparison of Dipole Parameters With and 
Without Extinction Template 



FIG. 3: Likelihood curves for different maximum redshifts 
in 2MRS. Any galaxies with \b\ < 8° are removed from the 
sample. 



resolution^ NSIDE=8 or 16. We adopt NSIDE=128 for 
the rest of this paper, except where otherwise noted, and 
also cut conservatively so that pixels with any area at all 
with \b\ < 8° are cut. 



C. Observational Constraints on Dipole Amplitude 
as a Function of Redshift 

With \b\ < 8° excised from the map (see Fig. [2|, we 
can apply the formalism outlined in Sec. |III| directly to 
different subsets of the survey. We pixelize the data us- 
ing HEALPix, meaning that we take the Galactic coor- 
dinates given in survey data and assign the given galaxy 
to the pixel corresponding to those coordinates. Fig. |3] 
gives the likelihood of the dipole amplitude, jC{A)^ for 
different subsets of the entire survey, reaching out to 
^max = 0.02,0.04,0.06,0.08, and 0.10 (where the last 
value in that list represents the entire survey except for 
25 sources at an assortment of higher redshifts). 

Note that the behavior is as expected in several re- 
gards: 

• The dipole amplitude A starts off larger and grows 
smaller as we go out further in redshift. 

• A converges to a certain value. This should happen 
simply because we run out of sources as we go to 
higher and higher redshift (e.g., 41,446 sources are 
at z < 0.06, while the total number of sources with 
z < 0.10 is 43,506.) 

• Because the redshifts in this sample are relatively 
low in cosmological terms, going out only to z ^ 
0.10, we expect that the dipole amplitude should 
remain on the order of 10~^, and it does. 



Resolution at some NSIDE roughly corresponds to pixel size ( 
l°x(60/NSIDE). 



If we wish to take Galactic extinction into account, the 
formalism we are employing to find dipole amplitudes 
and directions allows for very straightforward incorpora- 
tion of extinction due to dust as a systematic template. 
The 2MRS maps are already extinction-corrected, so this 
is much more of a sanity check than anything else. 

We use the SFD map [61], an intensity map of the 
sky at 100 microns and a reprocessed composite of the 
COBE/DIRBE and IRAS/ISSA maps. The resolution is 
on the order of a few arcminutes, and so the SFD map 
can in general be used to reliably derive extinction due 
to dust, assuming a standard reddening law. This works 
best away from the plane of the Galaxy, since within the 
Galactic plane dust conditions tend to fluctuate much 
more strongly on small scales than they do away from the 
Galactic plane (with the possibility of multiple dust tem- 
perature distributions, variable grain sizes, etc.). How- 
ever, since we work almost exclusively well away from the 
Galactic equator, we do not expect this to be an issue. 

The SFD map is nearly parity-even in Galactic coordi- 
nates, as the Galaxy is itself nearly parity-even, so when 
the sky is symmetrically cut, extinction is not likely to 
contribute to, or diminish, a dipole-like (parity-odd) sig- 
nal, at least not in the 2:-direction. However, we still in- 
clude it as a precaution against a known potential source 
of systematic error. 

The results are best presented in the form of a direct 
comparison, in Table HI Entries in this table take into 
account the entire 2MRS sample, 43,506 galaxies (before 
the symmetric cut in Galactic latitude) with 0.00 < z < 
0.10. 

The results change only slightly when the SFD tem- 
plate is added, and are statistically consistent with the 
no-dust results. The fact that the dipole amplitude drops 
slightly with addition of the template is an indication 
that a very small amount of the dipole power in the 2MRS 
map can be attributed to the pattern set up by the dis- 
tribution of dust in our Galaxy. 

The final row of Table |III| gives the results on dipole 
amplitude and direction when we include not only the 
SFD map but also the five £ = 2 and seven £ = 3 modes 
as templates. This ensures, as discussed earlier and elab- 
orated upon in Appendix |Bj that we are detecting only 
dipole signal and are not receiving contributions from 
higher multipoles. The contributions from higher multi- 
poles are almost negligible in this case, though they do 
slightly affect the direction of the detected dipole signal. 
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TABLE III: Comparison of dipole parameters with systematics templates vs. without templates, for 2MRS. The first column 
gives the Galactic latitude b of the cut; the second column identifies any systematic template present; the third gives the 
HEALPix NSIDE parameter; the fourth gives the number of sources that were still available after the cut was made; the fifth 
gives the dipole amplitude with highest likelihood; the sixth and seventh give / and b of the best-fit dipole; and the eighth 
and ninth give the 68 and 95 percent confidence intervals on the amplitude of the dipole. The first row gives results when we 
include no systematics templates; the second when we include only the SFD dust map as a template, and the third when we 
include SFD and quadrupole and octopole modes as templates. 



|6|> 


Systematics 


NSIDE 


N 


^peak 


I 


b 


68 percent CI 


95 percent CI 


8.0 


none 


128 


41834 


0.124 


228.0 


38.7 


0.116 - 0.132 


0.108 - 0.141 


8.0 


SFD Dust 


128 


41834 


0.118 


222.3 


38.3 


0.110 - 0.126 


0.102 - 0.135 


8.0 


SFD Dust + Quad + Oct 


128 


41834 


0.120 


213.8 


35.2 


0.111 - 0.128 


0.103 - 0.137 



E. Dipole Parameters as a Function of Sky Cut 

With the SFD extinction template in place, and again 
using the entire 2MRS sample out to z = 0.10, we may 
also vary the location where the cut in Galactic latitude 
is placed. Verifying that the placement of the cut (as 
long as it is kept at least as aggressive as the \b\ < 8° cut) 
does not affect the results beyond widening our error bars 
serves as a check for any source of systematic error that 
varies as a function of Galactic latitude. Most notably, 
any star-galaxy confusion that might creep into the sur- 
vey (very unlikely in the case of this particular survey) 
would vary strongly as a function of Galactic latitude, 
with the density of stars dropping precipitously as one 
moves away from the Galactic equator, and so this test 
serves to verify that star-galaxy confusion is not a major 
contributor to the detection of a dipole. It also helps to 
guard against the p ossibility that variations in sky cov- 
erage (see Sec. IV J) affect the dipole signal. (Sky cover- 
age is better at higher Galactic latitudes, since extended 
sources cannot be observed very close to very bright stars 
- less than 2 percent of the area of a typical high-latitude 
sky is masked by stars, as noted by Erdogdu et al. [17 .) 

Note that as the sky cut becomes more and more ag- 
gressive, we expect the error bars on the observed value 
of A to become wider (simply because we have less data 
and therefore looser constraints), but we expect the best- 
fit/peak-likelihood value of A itself to remain consistent 
with values found given less aggressive sky cuts. If A 
shifts in such a way that more aggressive sky cuts yield 
amplitudes inconsistent with amplitudes from maps with 
less aggressive sky cuts, this is a fairly good indication 
that star-galaxy confusion or another systematic effect 
that varies with Galactic latitude is at play. While the 
angular dimensions (and thus angular cuts) are not what 
determine the amplitude of the dipole, for very aggressive 
cuts we are left with far less data than we are with min- 
imal cuts, and this means that (a) the "measurement er- 
ror" on the observation becomes greater, and, even more 
importantly, (b) the cosmic variance associated with the 
theoretical prediction becomes much greater. Due to (b) 
in particular, a much wider range of peak-likelihood val- 
ues for the amplitude becomes consistent with theory, as 
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FIG. 4: Results for the dipole amplitude as a function of cut 
in Galactic latitude b. Here bcut indicates that \b\ < bcut was 
cut out of the map. Notice the very wide cosmic- variance 
band (yellow shaded region) around the theoretically pre- 
dicted value for the dipole amplitude, and the consistency of 
all observed values within cosmic- variance limits. The small 
green (blue) error bars indicate 68% (95%) measurement er- 
rors. The measurement errors are tiny in comparison with 
cosmic variance. Possible contamination from the quadrupole 
and octopole, which becomes more important as more of the 
sky is cut, is taken into account by including all ^ = 2 and 
i = 3 modes as systematics templates in the analysis. 



we cut the sky more aggressively. 

Results for different sky cuts are found in Fig. [4j The 
results are all consistent with one another, and with 
theory (all measurements being just outside the one- 
sigma cosmic- variance band) . Note that we have included 
quadrupole and octopole modes as systematics templates 
(see Appendix B) in creating this plot since £ = 2 and 
£ = 3 modes affect the error bars somewhat for more 
aggressive cuts. 



F. Dipole Amplitude: Theory vs. Observation 

Now that we have established basic consistency among 
dipole determinations with different sky cuts, we go back 
to the least restrictive cut, at \b\ < 8°, and keep the SFD 
template in place. We proceed to compare theory and 
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observation in dipole amplitude as a function of redshift 
in 2MRS. 

As noted earlier, it is important to take the bias of 
tracers of the LSS into account when producing theoreti- 
cal expectations for the clustering of these objects. Frith 
et al. [55] find the bias in the 2MASS Kg band to be 
1.39 ± 0.12, employing a technique that uses constraints 
on the galaxy power-spectrum normalization as well as 
previous constraints on ag. We therefore adopt 1.4 as the 
value of the bias for both 2MRS and 2MASS in general. 
The qualitative conclusions drawn from these surveys do 
not depend strongly on the precise value adopted for the 
bias. Note that, for constant bias 6, dipole amplitude is 
proportional to the bias, A (x b. 

There are cosmic- variance errors on all theoretical pre- 
dictions: 



SCi = 



(2^+l)/sky 



(26) 



Theoretical Dipole Amplitude 
Observed Dipole Amplitude 




0.10 



0.10 



We can relate Q to the amplitude A via Ci = (47r/9)A^ 
and so doing error propagation to get the cosmic- variance 
error on the amplitude, we have 



SA = 



1 



3 /sky 



-A. 



(27) 



This allows us to plot cosmic-variance uncertainties in 
both Ci and A. The basic result of doing so is shown 
in Fig. [5] Note that the errors on the observations are 
very small in comparison with cosmic variance. There- 
fore, for the rest of this section, we consider measurement 
errors negligible and focus only on how observational re- 
sults compare with theory within the bounds of cosmic 
variance. 

The dipole amplitude, both theory and measurement, 
decreases as redshift increases, exactly as it should given 
our previous arguments that averaging over more struc- 
ture at larger distances yields lower values of the dipole. 
Whether the amplitude is expressed as A or as Ci, the 
observational results are consistently lower than the the- 
oretical expectations. If these measurements for different 
^max were all independent, there would be a highly signif- 
icant inconsistency between theory and observation, but 
the measured values are highly correlated since samples 
with higher Zmax contain all samples with lower ^max- 

More specifically, correlations between a bin i going 
out to 2;max,i ^ud another bin j going out to Zmax,2 are 
calculated as 



Gov [Cii.Cjj] 



21 



(28) 



where we have made the usual assumption that the 
galaxy overdensity is a Gaussian random variable, so that 
Wick's theorem can be applied to obtain the expression 



FIG. 5: Top panel: Comparison of observations with theory 
for the dipole amplitude, as a function of how much of the 
2MRS sample is included (0.00 < z < Zmax)- All observed 
values are found including the SFD template, and with a cut 
at \h\ < 8°. For the purposes of calculating theoretical predic- 
tions, we take /sky = 0.86, corresponding to the |6| < 8° cut.) 
Bottom panel: The same results, only with the dipole power 
Ci rather than dipole amplitude A. In all cases, 68 percent 
error bars on observations are shown, but are either invisible 
or nearly invisible due to how tiny they are. See Sec. |III| for 
details of the procedure to convert between A and Ci . 



above. Then the correlation coefficient is 



Coy {Cjj.Cjj) 



VGov(a„a,)Gov(C,„C,,) 



(CiiCjj) 



(29) 



Gorrelations between 2MRS samples range from 0.42 (be- 
tween the full (0.00 < z < 0.10) sample and the smallest 
(0.00 < z < 0.01) sample), to 0.81 (between the full 
sample and the second-smallest (0.00 < z < 0.02) sam- 
ple), to well over 0.99 for many combinations of samples. 
(This is also the reason why all bins have similar signifi- 
cance as compared with one another.) Therefore, rather 
than being a 10-sigma inconsistency between theory and 
observation. Fig. [5] represents only slightly more than a 1- 
sigma discrepancy. The next section will find the precise 
"discrepancy" rigorously. 
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FIG. 6: Posterior probability of theoretical Ci given observed 
Ci as a function of the theoretical value, for the full sample of 
2MRS galaxies. The observed Ci determines the likelihood 
distribution for Cl^, and we can then compare the ACDM 
value (vertical magenta line) for Cl^ to that distribution. The 
ACDM value is clearly consistent with the observation. 



doing tomography in this way, and so we only perform 
this analysis on the full 2MRS sample, 0.00 < z < 0.10. 
The results are shown in Fig. [6j 

Whether we calculate a simple signal/noise ratio to 
compute the significance of Cl^ results given (7°^®, as 
in Fig. [5j or whether we use the more detailed compar- 
ison of Fig. [6j the qualitative conclusion is the same: 
the ACDM prediction matches observations within ap- 
propriate cosmic- variance limits. Note that the reason 
why these two strategies do not match up quantitatively 
is that the signal/noise strategy assumes cosmic variance 
is symmetric, while using the P(Cj^|C°^®) distribution 
takes into account the asymmetry of cosmic variance, 
particularly at the very low i at which we are working. 
This is also why significances in Ci and A do not match 
up with one another exactly. 



G. Comparison of Theory and Observation for 
Dipole Amplitude 

Given a Gaussian field on the celestial sphere with ob- 
served angular power spectrum the power is x^- 
distributed, and the likelihood of a given theoretical value 
is 



lnP(Cf|Cf^) = £ 



2e- 



^obs 



In 



^obs 



-InC; 



'obs 



(30) 

(see, for example, Chu et al. [62 )7 Here the observed 
quantity is treated as a realization of the theoretical 
value Cf^. For ^ = 1, this simplifies to 



\nP{Cl^\Cf'^ 



^obs C^^^ 



\nCf' (31) 



Again, we treat this as a likelihood, so that P is a func- 
tion of the theoretical model Cf^^ with the observed 
quantity held fixed. Then, as usual, we can plot the 
likelihood of a parameter value (in this case, theoreti- 
cal Cl^) and see where our "actual" theoretical Cl^ falls 
with respect to that distribution. In each redshift bin, we 
could generate a different likelihood distribution based 
on the observation, and then compare to the actual Cl^ 
in each case. However, because of the very high corre- 
lations between redshift samples, we gain very little by 



This expression can be derived by noting that the random vari- 

^obs 

able Y = (2i + 1)—^ is x^-distributed with 2^ + 1 degrees 

of freedom. Inserting this expression for Y into the general 
expression for a distribution, and then using the fact that 
P(Y)dY = P(Ci)dCi^ it is relatively straightforward to show 
that the proportionality for P{C() given in Chu et al. [62] holds, 
and from there the expression for the log-likelihood given above 
immediately follows. 



H. Observational Constraints on Dipole Direction 
as a Function of Redshift and Sky Cut 

Up to this point, we have been focusing on the dipole 
amplitude and comparing theoretical and observed am- 
plitudes. The direction of the dipole, however, is also a 
quantity of considerable interest. 

As discussed in Sec. [Hj there are three major types 
of dipoles that could contribute to any detected dipole 
in objects that trace large-scale structure: the local- 
structure dipole, the kinematic dipole, and the intrin- 
sic dipole. At the scales probed by 2MRS, we expect 
the local-structure dipole to completely dominate other 
contributions, since it is on the order of 10~^ while the 
kinematic dipole falls two orders of magnitude below this 
and the intrinsic dipole may very well fall even further 
below that. Therefore, there is no reason to expect that 
the direction of the 2MRS dipole should align with the 
direction of the CMB dipole, as we would expect it to do 
if the kinematic dipole were dominant at these scales. 

That said, it should not be at all surprising if the 2MRS 
local-structure dipole points somewhere near the CMB 
kinematic dipole. The reason for this has to do with 
what generates the motion that gives rise to the kine- 
matic dipole. As discussed in Sec.jlTj the total velocity of 
the Sun with respect to the Local Group is directed along 
almost the same line as the velocity of the Local Group 
with respect to the CMB rest frame, but in the opposite 
direction. So the direction of the Sun's total motion with 
respect to the CMB rest frame is essentially the same as 
the direction of the Local Group's motion with respect 
to the CMB rest frame, but the speed is lower than that 
of the Local Group since the contribution of the Sun's 
motion with respect to the Local Group gets subtracted 
off. The Local Group moves in a certain direction with 
respect to the CMB rest frame because of the gravita- 
tional pull of structure in the relatively nearby universe. 
The acceleration due to gravity of the Local Group, as 
determined via flux- weighted dipole measurements, is di- 
rected less than 20° away on the sky from the direction of 
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FIG. 7: Top panel: Likelihood associated with each dipole 
direction on the sky, marginahzed over amphtude, shown for 
2MRS redshift shells 0.00 < z < 0.01 (leftmost multicolored 
oval), 0.00 < z < 0.03 (uppermost multicolored oval), and 
0.00 < z < 0.10 (rightmost multicolored oval). We assume a 
\b\ < 8° cut, and incorporate the SFD dust systematic tem- 
plate and quadrupole and octopole templates. The color scale 
represents normalized likelihood as a function of direction. 
The single-colored disc that overlaps with one of the multi- 
colored likelihood ovals represents the direction of the CMB 
kinematic dipole, with error bars exaggerated to a circle of 2 
degrees in order to make the position clearly visible on the 
map. Bottom panel: Confidence intervals for the direction 
of the dipole in the full 2MRS survey, with the position of 
the CMB dipole shown. Agreement was not expected, but it 
is reassuring that the 2MRS projected dipole does lie in the 
same general region of sky as the CMB dipole. 



the velocity of the Local Group [29 . Therefore, insofar 
as the local-structure dipole gives information about the 
clustering of local structure and the direction of the accel- 
eration due to gravity of the Local Group, it is expected 
that it should point in at least the same general direc- 
tion as the CMB kinematic dipole, which is generated 
in part by the velocity of the Local Group induced by 
its acceleration due to gravity. Since the local-structure 
dipole is a 2D-projected quantity rather than one that 
preserves radial information, it is not a perfect indica- 
tor of where gravitational pulls on the Local Group are 
coming from. But we do expect the direction of the local- 
structure dipole to feel some influence from the direction 
of the CMB kinematic dipole. 

The observational results for the direction of the dipole 
are displayed in Fig. [7[ The results align with the quali- 
tative expectations detailed above. It turns out that the 



2MRS dipole direction is indeed not consistent with the 
direction of the CMB kinematic dipole, but still within 
the same basic region of sky. 

Note that the constraints on the dipole direction are 
actually tighter for the very small 0.00 < z < 0.01 sample 
than for the full sample of 2MRS galaxies, despite the 
fact that the number of sources is an order of magnitude 
smaller. This is not anomalous since the higher-redshift 
sources actually decrease the prominence of the dipole 
in local structure, producing a result with roughly half 
the total signal/noise as in the case where we take the 
^max = 0.01 subsample. 

The best-fit direction for the 2MRS dipole is (/,6) = 
(228.0°, 38.7°). Erdogdu et al. [17 find that the 2MRS 
number-weighted dipole (the quantity they analyze that 
is "closest" to our dipole) is at (/,6) = (218°, 33°) in the 
CMB frame, in close agreement with our results. 



I. Cutting the Supergalactic Plane 

As a final check, we wish to know how much of the 
dipole signal in 2MRS is coming from the vicinity of the 
supergalactic plane (SGP), a planar structure in the local 
galaxy distribution (Lahav et al. [63 ). We therefore pro- 
gressively excise more and more of the SGP and see how 
much the amplitude of the dipole dies away. We com- 
pare this to the effect of excising similar areas from the 
vicinity of the supergalactic poles. We expect that there 
should be more sources near the supergalactic plane, and 
that the dipole should die away much more quickly when 
the supergalactic plane is excised than when similar areas 
around the supergalactic poles are excised. 

This check will become more important as we pro- 
ceed to perform our analysis on surveys that probe much 
larger radial distances than does 2MRS, as the structure 
associated with the supergalactic plane will only con- 
tribute to the dipole on relatively nearby scales, and the 
effect should diminish as we probe to larger and larger 
redshifts. 

We find the results in Table HVl when we use the usual 
\b\ < 8° Galactic cut and also include a cut around the su- 
pergalactic plane (|SGB| > in the first column^) or a cut 
around the supergalactic poles (|SGB| < in the seventh 
column). The table is aligned so that cuts around the 
supergalactic plane are in the same row as cuts of similar 
area around the supergalactic poles. The /sky columns 
give the fraction of the sky that remains when we per- 
form the given cut. We also calculate /sources ^ the fraction 
of the total number N of sources that remain when we 
perform the given cut. The cuts in SGB (less than 2.0, 



Supergalactic coordinates SGB and SGL are defined in analogy 
to Galactic b and /, where SGL is the azimuthal coordinate, and 
SGB=0 corresponds to the middle of the supergalactic plane, so 
we care only about making cuts in SGB. 
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TABLE IV: Key patterns in cutting in supergalactic coordinates, for 2MRS. 
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/sky 




/sources 
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J sky 
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/sky 
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0.86 


41834 


1.00 
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2.0 


0.82 


39964 


0.96 


1.16 


0.12 


74.82 


0.84 


41234 


0.99 


1.17 


0.12 


5.0 


0.78 


37124 


0.89 


1.14 


0.12 


65.90 


0.81 


39867 


0.95 


1.18 


0.11 


10.0 


0.70 


32673 


0.78 


1.11 


0.12 


55.73 


0.74 


36882 


0.88 


1.20 


0.09 


20.0 


0.55 


24799 


0.59 


1.08 


0.11 


41.15 


0.59 


30321 


0.72 


1.22 


0.05 



All Sources, Equatorial Coordinates 



Ks<13.5, Equatorial Coordinates 




All Sources, Galactic Coordinates 



Ks<13.5, Galactic Coordinates 




FIG. 8: Left side: All 2MASS sources, in equatorial and Galactic coordinates. Note the very strong discontinuity in the selection 
function, visible in both images (especially in the top image, where it appears as a horizontal line), at declination around 20 
degrees. Right side: 2MASS sources with Ks-hand magnitude less than 13.5. The survey has nearly uniform completeness 
when this criterion is imposed. 



5.0, 10.0, 20.0 degrees; greater than 74.82, 65.90, 55.73, 
41.15 degrees) were chosen so that equal areas around the 
plane and around the poles would be cut if there were no 
cut in G alact ic b (as in our tests on the BATSE cata- 
log; Sec. |VII| ). Since there is a cut in Galactic latitude 
here, /sky does not match up exactly between the cuts 
around the supergalactic plane and poles, but the values 
are still close, and in any case they are normalized in the 
/sources //sky Calculation; we employ the same cuts for all 
surveys tested in this paper. 



From Table IV, we learn the following: 



• The fraction of sources associated with the super- 
galactic plane is greater than the fraction of sources 
associated with the supergalactic poles for every 
/sky. This is our first indication that a greater- 
than-random portion of the dipole signal comes 
from the vicinity of the supergalactic plane. 

• The ratio /sources //sky, which gives a measure of 
how overdense the uncut portion of the sky is, 
dwindles steadily as we cut more and more of the 
SGP, but increases as we cut more and more area 
around the supergalactic poles. Therefore, as ex- 
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pected, there are more sources near the supergalac- 
tic plane than near the supergalactic poles, and this 
is true essentially regardless of how much of the 
area around the plane/poles is cut. 

In general, we conclude that more of the dipole signal 
comes from the area of the supergalactic plane than from 
the vicinity of the supergalactic poles. This serves as a 
good check that the source of the dipole signal in the 
relatively local structure surveyed by 2MRS is generally 
where we expect it to be. When we perform analyses of 
higher-redshift objects in Sec. |VII| and |VIII[ the super- 
galactic plane should not "show up" as it has here. 



J. 2MRS: Conclusion 

We conclude our analysis of the dipole signal in 2MRS 
by pointing out that all results are consonant with theo- 
retical expectations. This comes as no surprise given that 
2MRS was the most well-controlled and well-understood 
of the surveys we analyze here, and was being treated 
exhaustively as something of a model for our other anal- 
yses. 

It should also be noted that we have verified that the 
results above do not change appreciably when corrections 
are made for the 2MASS coverage map (see Sec. [V|). 

We now proceed to apply our dipole analysis to the 
fuh 2MASS dataset. 



V. DIPOLE IN 2MASS 

We now analyze the full 2MASS survey in a manner 
similar to how we analyzed the very well-characterized 
2MRS subsample. The challenges associated with ana- 
lyzing 2MASS as a whole are greater, in part because 
the full sample of 2MASS galaxies (1.6 million extended 
sources) does not have uniform completeness across the 
entire sampled sky. 

It should be noted that some previous results concern- 
ing the 2MASS galaxy distribution stand in some ten- 
sion with ACDM predictions. For example. Frith et al. 
[BT point out that the angular correlation function and 
angular power spectrum of 2MASS galaxies (under cuts 
reasonably similar to the ones we perform here) display 
fiuctuations that are 3-5 sigma out of line with ACDM 
predictions. We focus attention here on the dipole alone, 
which of course sacrifices a certain amount of informa- 
tion with respect to what could be gained from analysis 
of the entire power spectrum, but also lends itself to much 
better and more detailed analysis of contributions to the 
signal at this one multipole. 



A. Selection Cuts 

We make several cuts to the sample of 2MASS galaxies 
in order to ensure uniformity of the sample: 

1. As shown in Fig. [8) the biggest issue in connec- 
tion with survey completeness is that the selection 
function has a sharp discontinuity for galaxies with 
i^5-band magnitude greater than roughly 13.5. We 
therefore cut out all these sources, roughly 2/3 of 
the sample, at the outset, and consider only that 
portion of the survey with nearly uniform complete- 
ness over the entire sky (with the exception of the 
Galactic plane), that is, sources with Kg < 13.5. 

2. We must make a more aggressive Galactic cut than 
we did for 2MRS in order to ensure that star-galaxy 
confusion does not come into play. Mailer et al. 
[65] and Skrutskie et al. [58] note that the 2MASS 
XSC is highly reliable and complete for \b\ > 20° 
(more than 98 percent galaxies rather than stars 
at these latitudes), but that star-galaxy confusion 
is an increasingly large problem at lower latitudes: 
the XSC is 10 percent stars for 5° < \b\ < 20°; and 
within the Galactic plane, \b\ < 5°, there is addi- 
tional contamination by artifacts (10 to 20 percent) 
and Galactic extended sources (~ 40 percent) in- 
cluding globular clusters, open clusters, planetary 
nebulae, and giant molecular clouds (Jarrett et al. 
[66 ). In particular. Mailer et al. cross-correlate the 
2MASS stellar density ngtar with the XSC galaxy 
density as a function of the latitude of a symmetric 
(in Galactic coordinates) cut and find that includ- 
ing XSC objects with \b\ < 15° gives a galaxy-star 
cross-correlation that is higher in amplitude than 
the galaxy-galaxy autocorrelation, suggesting the 
presence of multiple-star systems mistakenly iden- 
tified as extended sources. However, this excess 
signal goes away for a cut of \b\ < 20°. Cutting at 
\b\ < 20° ensures less than 2 percent contamination 
from Galactic sources [55l |65] . 

3. We use the XSC confusion fiag (cc_f Ig) to elim- 
inate known artifacts (diffraction spikes, meteor 
streaks, infrared airglow, etc.). 

4. Again following Frith et al. [55], as well as Mailer 
et al. [65 , we also cut out bright {Ks < 12.0) ob- 
jects with (J — Ks) colors that are outside the range 
[0.7, 1.4] (see Bilicki [36 for alternative choices of 
which ranges to cut). This is a conservative mea- 
sure designed to get rid of a final set of objects 
which are in the 2MASS XSC but which are not 
extragalactic sources. This removes a few thousand 
sources. 

5. As explained later, we take the 2MASS sky cov- 
erage into account. The XSC does not have com- 
pletely uniform sky coverage given the presence of 
bright stars and other foreground objects that make 
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FIG. 9: Plot of the radial distribution of 2MASS galaxies as a 
function of redshift. Three different mean redshifts are shown, 
one which corresponds to the 2MRS distribution (the actual 
2MRS dN/dz is plotted in blue), and the other two of which 
correspond to photometric cuts in the full 2MASS survey of 
Ks < 12.5 and Ks < 13.5. 



it more difficult for the telescopes to detect ex- 
tended sources in particular directions. Although 
the pattern of sky coverage is parity-even (follow- 
ing the shape of the Galaxy) and unlikely to mimic 
a dipole in any way (as implied by the 2MRS anal- 
ysis), we still take this into account in the present 
analysis. 

6. K-corrections (corrections to magnitudes in a given 
passband that are made necessary by the fact that 
light can redshift into or out of a given range of 
wavelengths) for the Kg band can make a non- 
negligible difference in the calculation of a flux- 
weighted dipole or other quantity that depends on 
specifics of photometry. However, in this case we do 
not need to take them into account, because they 
are actually accounted for in our predictions: K- 
corrections are tied to the same (pseudo-)Doppler 
effect that helps to generate the kinematic dipole 
(see Sec.|ll|), and so accounting for them in obser- 
vational results as well would amount to double- 
counting. 

All photometric cuts are applied to 2MASS isophotal 
magnitudes - not total magnitudes, which are an extrap- 
olated quantity and viewed as less reliable for the pur- 
poses of this kind of analysis. While many analyses which 
use 2MASS data actually use the extrapolated magni- 
tudes (since, according to Jarrett et al. [67 , the isophotal 
magnitudes underestimate total luminosity by 10 percent 
for early- type and 20 percent for late- type galaxies), we 
stick here with the more conservative isophotal magni- 
tudes, especially since the cut at i^g < 13.5 is much surer 
to accomplish its purpose if the more conservative mag- 
nitude estimates are used. It is worth noting that the 
2MRS team used isophotal magnitudes in their sample 
selection (Huchra et al. [60^). 



B. Radial Selection Function 

With these photometric cuts applied, we are ready to 
proceed with the analysis. While no spectroscopic red- 
shifts are available for the 2MASS XSC as a whole, and 
the photometric redshifts that do exist are not particu- 
larly reliable, considerable information is available about 
the overall radial distribution of 2MASS galaxies. In par- 
ticular. Frith et al. [55 and others give the 2MASS radial 
selection function as 

dN ^ , 3^2 / /1.412^\^/^\ , , 

with z = 0.074 for Ks < 13.5 and z = 0.050 for 
Ks < 12.5. With these values of we can determine the- 
oretical predictions for the local-structure dipole (which 
is still dominant by two orders of magnitude over other 
contributions to the dipole at these scales) for these two 
photometric cuts. Combined with the 2MRS sample, 
which follows this same form for the selection function 
quite closely (see Fig. |9| and corresponds to approxi- 
mately z = 0.028, we can perform a comparison of the- 
ory and observation for multiple subsamples of the entire 
2MASS catalog. 

For Ks < 12.5 without (with) the \b\ < 20° cut, there 
are 127,030 (89,980) galaxies. For Ks < 13.5 without 
(with) the cut, there are 542,201 (381,586) galaxies. 



C. Systematic Checks: Extinction 

In 2MRS, extinction corrections were already applied 
to the magnitudes of the galaxies, but in 2MASS, the 
catalog values for the magnitudes are not corrected for 
extinction. This means that it becomes much more im- 
portant in this case to make sure that we have adequately 
controlled for the effects of extinction. Knowing the mag- 
nitudes is important to determine which objects get into 
the sample in the ffist place. We find that several thou- 
sand galaxies that do not make the Ks < 13.5 cut before 
extinction correction do make the cut when magnitudes 
are corrected for extinction.^ 

We have performed various extinction corrections, ex- 
perimenting with slightly different extinction coefficients 
R = Av/E{B - V) for the 2MASS Ks band (0.367 from 
Ho et al. [11]; 0.302 from the analogous UKIRT value in 
Schlafly and Finkbeiner [68j; cf. 0.35, which is used by 
Erdogdu et al. [17 , following Cardelli et al. [69]). We find 
that the results when extinction corrections are applied 
directly are essentially identical to the results obtained 



Note that extinction corrections always bring magnitudes down, 
since sources appear dimmer due to extinction, and so are as- 
signed higher brightness/lower magnitude when corrected for ex- 
tinction. 
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TABLE V: Comparison of dipole parameters without any templates, with SFD template, and with SFD, quadrupole and 
octopole templates, for 2MASS, for two different limiting iT-band magnitudes. 



For objects with Ks < 13.5: 


1^1 > 


Template 


NSIDE 


N 


-^peak 


/ 


b 


68 percent CI 


95 percent CI 


20.0 


none 


128 


386008 


0.089 


303.4 


7.3 


0.086 - 0.092 


0.083 - 0.095 


20.0 


SFD 


128 


386008 


0.088 


305.0 


4.5 


0.085 - 0.091 


0.082 - 0.094 


20.0 


SFD + Quad + Oct 


128 


386008 


0.104 


268.4 


0.0 


0.100 - 0.108 


0.096 - 0.112 


For objects with Kg < 12.5: 


1^1 > 


Template 


NSIDE 


N 


-^peak 


/ 


b 


68 percent CI 


95 percent CI 


20.0 


none 


128 


91008 


0.0848 


275.0 


28.2 


0.078 - 0.091 


0.072 - 0.097 


20.0 


SFD 


128 


91008 


0.0812 


276.3 


25.9 


0.075 - 0.088 


0.069 - 0.094 


20.0 


SFD + Quad + Oct 


128 


91008 


0.134 


267.3 


8.5 


0.126 - 0.142 


0.117 - 0.150 



when the SFD dust systematic template is applied as it 
was to the 2MRS maps, so we explicitly present only the 
SFD-template results here. 

Table [Vl shows the effects of the inclusion of the SFD 
template for the Ks < 13.5 and Ks < 12.5 maps. As 
should be clearly visible from a quick glance at the values 
in the table, very little changes when the SFD template is 
included, and there is substantial overlap even of the 68 
percent confidence intervals for each of the no-template 
cases with each of the corresponding SFD-template cases. 
We conclude that although the results shift slightly, and 
therefore it is worth keeping the SFD template in our 
analysis, extinction does not have a substantial impact 
on the dipole results. This is as expected based on con- 
siderations of how extinction affects 2MASS coverage, as 
outlined in Jarrett et al. [70 and Jarrett et al. [66 : the 
completeness of the 2MASS XSC is much more adversely 
affected by source confusion than by extinction. 

In the same table, we also show the effect of includ- 
ing templates for the quadrupole and octopole modes, 
and note that including these templates shifts the peak 
dipole amplitude by more than 10 percent for Ks < 13.5 
and more than 50 percent for Ks < 12.5. (Note also the 
changes in direction in each case, which actually serve 
to bring into agreement the dipole directions for the two 
different photometric cuts.) Based on the results pre- 
sented in Appendix [Bj it is not surprising that inclusion 
of these templates has more of an effect than it did for 
2MRS, since our sky cut is significantly more aggressive 
here. 



D. Systematic Checks: Coverage Map 



makes observation of distant background galaxies impos- 
sible. Sky coverage, which ranges from to 1 within 
a given pixel, tends to be well above 0.98 for the high- 
Galactic-latitude sky. Data products from 2MASS in- 
clude coverage maps^^ that indicate coverage as a func- 
tion of direction. 

We convert these maps into the same HEALPix pix- 
elization scheme we use to pixelize all the surveys in this 
paper, including 2MASS itself. Each HEALPix pixel con- 
tains at least 4, and up to 19, "subpixels" associated with 
the pixelization of the 2MASS coverage maps, so resolu- 
tion is not an issue. See Fig. [lO] 

There are several ways in which we could take these 
coverage maps into account in our analysis. First, we 
could mask out all HEALPix pixels that have an average 
coverage less than some threshold (the threshold is usu- 
ally chosen as 0.98 in the literature; see, e.g.. Ho et al. 
[31]). Second, we could mask out all pixels that have any 
subpixel with coverage less than some threshold. Third, 
we could use the entire coverage map as a systematic 
template. We have not found a case in which it makes 
anything even close to a statistically significant difference 
which of these strategies we choose, so we choose the op- 
tion that is simplest and arguably best: we use the entire 
coverage map as a systematic template. This has the ad- 
vantage of not privileging any particular threshold, but 
rather taking the variation in coverage over the entire 
sky into account evenhandedly. This accounts for the ac- 
tual pattern of observations on the sky and weights them 
accordingly (cf. our treatment of the BATSE exposure 
function in Sec. |VII[ ). In any case, given that the coverage 
map closely follows the shape of the Galaxy (foreground 



Uniform completeness in an infrared survey like 

2MASS is impossible due to the presence of foreground lO See http : //www ipac . caltech . edu/2mass/releases/allik^ 

stars. In some directions, the presence of foreground stars doc/sec2_6f .htrnTI 
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Average Coverage in Each Pixel, Galactic Coordinates 



o.s 




FIG. 10: 2MASS sky coverage as a function of direction for 
the entire sky. Coverage can be zero (full masking) or one (no 
masking), or anywhere between these extremes. We use this 
map as a systematic template in the dipole formalism. 



stars are, after all, our primary concern), and the Galaxy 
is nearly parity-even, we do not expect coverage to con- 
tribute significantly to the (parity-odd) dipole anyway; 
results are very much in accord with these expectations. 

In summary: We apply the coverage map throughout 
this section (as we did, in fact, for our 2MRS analysis as 
well), but find that the results do not change appreciably 
as a result. 



E. Systematic Checks: Sky Cut and Supergalactic 
Plane 



We perform the same cuts in Galactic latitude as we 
did for 2MRS, but with the expectation that the most 
reliable results will come for the \b\ < 20° cut rather 
than \b\ < 8° as in the case of 2MRS. See Fig. fl2]for a 



visual capture of the results, with cosmic variance on the 
theoretical prediction taken into account. 

The basic conclusion here, as in the case of 2MRS, is 
that in no case are results outside of the limits expected 
given cosmic variance. There is tension between the 
^cut = 10° sample and other samples for both Kg < 13.5 
and Ks < 12.5, but that is to be expected given the much 
higher potential for star-galaxy confusion when such low 
latitudes are included. The same could be said of the 
^cut = 15° sample for i^g < 13.5. In any case, cosmic 
variance dominates the error budget for all cases, and 
dependence of results on Galactic cut indicates no seri- 
ous contamination from star-galaxy confusion or other 
systematic effects that vary with Galactic latitude for 
samples with Galactic \b\ < 20°, so we follow Skrutskie 
et al. [58 , Frith et al. and others in taking this as 
our fiducial cut. 

We also perform the same test as in the 2MRS case 
where we cut in supergalactic latitude SGB. We create 
summary tables giving the key patterns, as we did for 
2MRS in Table IV Results are presented for Kg < 13.5 
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FIG. 11: Top panel: Results for the dipole amplitude in the 
2MASS survey, as a function of mean redshift z of the galaxy 
sample. The z = 0.028 sample corresponds to the 2MRS 
galaxies with \b\ < 8° cut out; the z = 0.050 sample to Kg < 
12.5 with \b\ < 20° cut out in the 2MASS XSC; the z = 
0.074 sample to Kg < 13.5 with \b\ < 20° again cut out in 
the 2MASS XSC. Cosmic variance (yellow band around the 
theoretical prediction) is shown in the most pessimistic \b\ < 
20° case for all of the samples. Bottom panel: Comparison 
of the ACDM value for Cl^ with the expected distribution 
of Cl^ given the observed value for the Ks < 12.5 sample. 
Both panels demonstrate sound agreement between theory 
and observation for all three subsamples of 2MASS. 



and Ks < 12.5 in Table |VI[ The observations that we 
made for 2MRS again (generally) hold here: in partic- 
ular, /sources //sky dccrcascs mouotouically as more and 
more of the supergalactic plane is excised, and the ra- 
tio increases (almost) monotonically as more and more 
of the area around the supergalactic poles is excised. 



F. Dipole Amplitude as a Function of 
Redshift/Photometric Cuts 

Taking as most reliable the case with a cut for \b\ < 
20°, and keeping the SFD, quadrupole, and octopole tem- 
plates in place, we proceed to compare theoretical pre- 
dictions with observational results for the dipole ampli- 
tude for the two different magnitude cuts we have used, 
Ks < 13.5 and Ks < 12.5, which correspond to ^ = 0.074 
and 0.050, respectively (where, recall, z is defined un- 
der Eq. ([32|)). Again, the 2MRS sample corresponds to 
z = 0.028, and we include this data point in our compar- 
isons as well. 



See Fig. 11 for results. Note that measurement errors 
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TABLE VI: Comparison of dipole parameters when performing various cuts in supergalactic coordinates, for 2MASS, for two 
different limiting K-hand magnitudes. 



For objects with Kg < 13.5: 


|SGB| > 


/sky 


N 


/sources 


/sources 
/sky 


^peak 


|SGB| < 


/sky 


N 


/sources 


/sources 
/sky 


^peak 


0.0 


0.65 


386008 


1.00 


1.53 


0.10 














2.0 


0.63 


368077 


0.95 


1.52 


0.09 


74.82 


0.65 


385740 


1.00 


1.53 


0.10 


5.0 


0.59 


342390 


0.89 


1.51 


0.08 


65.90 


0.64 


378368 


0.98 


1.53 


0.10 


10.0 


0.52 


301028 


0.78 


1.50 


0.08 


55.73 


0.60 


355464 


0.92 


1.53 


0.07 


20.0 


0.39 


225502 


0.58 


1.49 


0.13 


41.15 


0.49 


295090 


0.76 


1.55 


0.07 


For objects with Kg < 12.5: 


|SGB| > 


/sky 


N 


/sources 


/sources 
/sky 


^peak 


|SGB| < 


/sky 


N 


/sources 


/sources 
/sky 


^peak 


0.0 


0.65 


91008 


1.00 


1.53 


0.13 














2.0 


0.63 


86657 


0.95 


1.52 


0.11 


74.82 


0.65 


90951 


1.00 


1.53 


0.13 


5.0 


0.59 


80177 


0.88 


1.50 


0.10 


65.90 


0.64 


89203 


0.98 


1.53 


0.12 


10.0 


0.52 


70212 


0.77 


1.48 


0.10 


55.73 


0.60 


83805 


0.92 


1.53 


0.10 


20.0 


0.39 


52528 


0.58 


1.48 


0.14 


41.15 


0.49 


69923 


0.77 


1.55 


0.10 



are once again tiny in comparison with cosmic-variance 
errors. 

The magnitude of the dipole in 2MASS has, of course, 
been calculated previously as part of computations of 
the entire power spectrum. In particular. Frith et al. [55] 
give Ci ^ 0.004 for both Ks < 13.5 and Ks < 12.5. 
This value converts to A ^ 0.054; cf. our values of 
A = 0.104 ± 0.004 and A = 0.134 ± 0.008 (68 percent 
confidence) for Ks < 13.5 and Ks < 12.5. While both 
our measurement and the Frith et al. measurement are 
in two-sigma cosmic- variance agreement with theory, the 
mutual discrepancy between these two observational re- 
sults is noteworthy. While we do not fully understand 
this discrepancy with the Frith et al. result, we note that 
our results are in better agreement if we do not marginal- 
ize over the dust map and the quadrupole and octopole 
templates. 



G. Dipole Direction as a Function of 
Redshift/Photometric Cuts 

In Fig. [13) we present the results for the dipole di- 
rection in both the Ks < 13.5 and more conservative 
Ks < 12.5 cases. Once again, dipole amplitudes are on 
the order 10~^ while the kinematic dipole is expected to 
be on the order 10~^, so no particular agreement with 
the direction of the CMB dipole is expected. We do ex- 
pect the Ks < 12.5 sample to give a direction relatively 
close to that of the 2MRS dipole, given the overlap in 
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FIG. 12: Top panel: Dipole amplitude as a function of bcut 
for 2MASS sources with Kg < 13.5. Bottom panel: Same 
for 2MASS sources with Kg < 12.5. See Fig. |t] for a fuller 
discussion of the significance of this type of plot. 
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FIG. 13: Confidence intervals to go witii results for the dipole 
direction in the 2MASS XSC, Kg < 13.5 (smaller circles) and 
Ks < 12.5 (larger circles). In both cases we apply a cut elim- 
inating \b\ < 20°, and apply the SFD map and quadrupole 
and octopole maps as systematic templates. The CMB kine- 
matic dipole direction is indicated. Like 2MRS, 2MASS is too 
shallow to expect agreement between its dipole direction and 
the direction of the CMB dipole, so this is not an anomalous 
result. 



the samples, and Kg < 13.5 to give a result close to 
Kg < 12.5. In fact, we would regard it as anomalous 
if the results were not all consistent with one another, 
if the samples were genuinely sampling the same popu- 
lation - larger samples would simply have smaller error 
bars than smaller samples. However, the populations be- 
ing sampled are different, given that the structure asso- 
ciated with z = 0.028, 0.050, and 0.074 are quite distinct 
from one another. So even internal inconsistencies (be- 
tween different values of z) are tolerable. Indeed, we find 
that the direction of the Kg < 13.5 dipole is not fully con- 
sistent with the direction of the Kg < 12.5 dipole, but 
especially when we include the quadrupole and octopole 
modes as systematics templates, as we do in Fig. 
inconsistency is very mild. 
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H. 2MASS: Conclusion 

We draw the basic conclusion that there are no anoma- 
lous results in applying tests of dipole amplitude and di- 
rection to subsets of the 2MASS dataset and comparing 
these results with theoretical predictions. We now pro- 
ceed to more critical tests using higher-redshift objects 
that might begin to probe the kinematic dipole, as these 
objects exist at scales on which the local-structure dipole 
should have become comparable to, or even smaller than, 
the kinematic dipole. 



VI. DIFFICULTIES WITH USING X-RAY 
SURVEYS 

We begin by considering surveys that detect very high- 
energy photons (X-ray and gamma-ray), and then ad- 
dress the opposite end of the spectrum (radio). 



First, a brief note on X-ray data. Flux- weighted 
dipoles have been previously calculated using the soft X- 
ray band (< 2 keV) data from ROSAT and the hard X- 
ray (2-10 keV) background data as observed by HEAOl- 
A2. For example, Plionis and Georgantopoulos [71 use 
the 1.5 keV (- 0.8 nm) ROSAT All-Sky Survey (RASS) 
data to calculate a flux-weighted dipole, and Scharf et al. 
[72 perform a similar analysis on HEA01-A2. Contami- 
nation is a major issue for both analyses. At the time of 
publication of these studies, only roughly three-quarters 
of the unresolved X-ray flux in the soft and hard bands 
had been accounted for (by extrapolation of objects re- 
solved in deep fields) [72], and theoretical modeling of 
populations contributing to the unresolved flux (AGNs, 
starburst galaxies, hot IGM in rich clusters, etc.) re- 
mained difficult. In the case of hard X-rays, Scharf et 
al. argue that at least a third of the structure in the data 
may be Galactic in origin (associated especially with the 
bulge), and soft X-rays are even more strongly contami- 
nated by Galactic emission. Treyer et al. [73] point out 
that in the soft band. Galactic emission is present as a 
contaminant at all scales; the hard band is better, but 
Treyer et al. still rely on the Galactic hard-band emis- 
sion model of Iwan et al. [7? to predict that for Galactic 
latitude > 20°, the variations in flux due to Galactic 
emission are less than 3 percent. 

All told, it is very difficult to remove the foreground 
in X-ray all-sky surveys successfully without some rel- 
atively uncertain modeling, without fairly serious suspi- 
cion of contamination, and without removing a good deal 
of the background too (especially in the soft band; Plio- 
nis and Georgantopoulos [71 estimate that Virgo con- 
tributes as much as 20 percent to the dipole amplitude 
in RASS). Theoretical predictions would also be difficult 
to make without a well-understood redshift distribution, 
especially given that the populations contributing to the 
X-ray background are not especially well-modeled. Hence 
we stay away from attempting to perform X-ray analyses 
here. 

Some authors have searched for a dipole in the XRB 
to the extent that it is possible to do so. The results that 
Scharf et al. find in hard X-rays for the flux-weighted 
dipole basically align with theoretical predictions of what 
they refer to as the Compton- Getting effect (another 
name for the kinematic dipole, following a paper by 
Compton and Getting [75 on the effect on cosmic-ray in- 
tensity of Earth's motion through the Milky Way). More 
recently, Boughn et al. [76j analyzed the same (HEAOl- 
A2) dataset and found a limit (95 percent confidence) 
on the amplitude of any intrinsic dipole at 5 x 10~^. 
However, given the difficulty of definitively separating 
extragalactic from foreground/Galactic emission in this 
dataset, and other problems already noted, significant 
uncertainty attends any analysis in X-rays, so we note 
all these results without making heavy use of them in 
the remainder of this paper. Similarly, other populations 
of objects detected at the very high-energy end of the 
spectrum, including blazars and clusters of galaxies de- 
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tected with gamma-ray satellites (see, e.g., Ando et al. 
[77 ), may be good targets for dipole searches in the long 
term, especially once their bias is better-understood and 
future surveys provide better statistics for the given tar- 
get population; but we do not pursue those here. 



VII. DIPOLE IN BATSE GAMMA-RAY BURSTS 

For analysis of a 2D-projected dipole (i.e., not flux- 
weighted; see Sec. |ll|) in very high-energy surveys, we 
turn instead to gamma-ray bursts (GRBs). GRBs are the 
most powerful explosions known in the universe, though 
their exact nature and progenitor objects remain under 
some debate, and their redshifts are difficult to mea- 
sure since GRB observations are not well-localized (see 
Fig. 14) and redshifts can only be measured from their 



afterglows, which must be matched up with the position 
of the original GRB, a highly nontrivial task given the 
error bars on the position of a typical GRB. A review 
of previous research on the dipole in the GRB distribu- 
tion, as well as presentation of results using the formalism 
outlined in Sec. |III|an d applied to the BATSE catalog, is 
presented in Sec. |VII[ In Sec. |VIII[ we move to the low- 
frequency end of the spectrum and present results from 
the NRAO VLA Sky Survey (NVSS). 



A. Previous Work on the Isotropy of GRBs 

Up through the mid-1990s, there was a long history 
of assessing the isotropy of gamma-ray bursts (GRBs) 
in an attempt to infer whether they were cosmological 
or Galactic sources (or a combination thereof). For ex- 
ample, Maoz [78] argued in 1993 that gamma-ray bursts 
could be shown to exist in an extended Galactic halo, 
some 130-270 kpc away from Earth, by detecting slight 
but well-defined deviations from spherical symmetry pre- 
dicted for such a halo population. While his analysis 
did indeed suggest that GRBs were nearby intergalactic 
objects, he argued that comparison with more specific 
models would be necessary before considering the case 
closed. 

In a similar spirit, Briggs et al. [79] argued persua- 
sively that the population of GRBs could not be Galac- 
tic, based on their observed isotropy. This study found 
that the Galactic dipole and quadrupole moments (calcu- 
lated very straightforwardly as (cosO) and (sin^ b — 1/3)) 
did not differ significantly from those predicted for an 
isotropic distribution. The majority of GRB models that 
assumed GRBs are a Galactic population were found 
to be in > 2a tension with the detected dipole and 
quadrupole moments, and hence the conclusion of this re- 
search was that GRBs are more isotropic than observed 
Galactic populations, suggesting either a nearby inter- 
galactic or, more likely, cosmological source. 

Scharf et al. ^] computed a fluence-weighted dipole 
(where fluence is flux integrated over the timespan of 



the burst) in analogy to the fiux-weighted dipoles dis- 
cussed in previous sections of this paper. Combining 
fluence-weighted dipole information with straightforward 
2D-projected dipole measurements (i.e., including pho- 
ton count information) better distinguishes a velocity 
dipole (due, as usual, to the Doppler effect and relativis- 
tic aberration) from other possible sources of anisotropy. 
This kind of test can be regarded as a supplement to the 
kinds of tests we perform here. 

The current consensus that GRBs are cosmological is 
based not only on the considerations discussed above and 
the absence of even a weak band corresponding to the 
Milky Way in the GRB distribution (Tegmark et al. [81j; 
see also Balazs et al. [82 ), but also (and especially) on the 
observation of optical. X-ray, and radio counterparts to 
GRBs that are clearly extragalactic (e.g., [83-86 ). Given 
the extragalactic origins of GRBs, we should expect a 
GRB dipole sourced by the same effects that give rise to 
a dipole in other sources we have analyzed. As far back 
as the mid-1990s, Maoz [87] predicted the dipole in the 
clustering of GRBs, combining the effects of relativistic 
aberration and the Doppler effect. These estimates are 
somewhat uncertain, but we too are unsure how precisely 
GRBs trace large-scale structure (e.g., if they can be de- 
scribed as a single population with a single bias, and so 
forth). Maoz finds that the amplitude of this (kinematic) 
dipole is of order A ~ 10~^ (to within uncertainties of a 
factor of two), which is still an order of magnitude larger 
than the CMB dipole, but much closer than the sources 
we dealt with in the previous section. Maoz estimated 
that a large (of order 10^) sample of GRBs would be 
necessary to detect the predicted dipole, and given that 
current catalogs offer only on the order of 10^ bursts, we 
do not expect an unequivocal detection. That said, we 
nevertheless run our tests, given that useful constraints 
can be placed on the maximum possible dipole amplitude 
even if we cannot confidently detect a dipole in currently 
available GRB catalogs. 



B. BATSE Data and Systematic Effects 

We perform our tests on GRBs in the BATSE catalog, 
data taken by the Compton Gamma-Ray Observatory. 
The BATSE (Burst and Transient Source Experiment) 
instrument onboard CGRO detected gamma-ray bursts 
within the nominal range of 50 to 300 keV. Other GRB 
datasets are available, including those from SWIFT and 
Fermi (known previously as GLAST), but we use the 
BATSE catalog because it has the most sources for an 
all- sky survey. 

The most recent catalog of GRBs from BATSE is the 
"current" catalog with 2702 sources. In this catalog, 
there are several complicating details that are worth not- 
ing in regard to what bursts make it into the catalog and 
which do not. 

BATSE employed scintillators sensitive to gamma rays 
from ~ 25 to 2000 keV. A burst triggered the instrument 
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BATSE GRBs m Galactic Coordinates 




BATSE Exposure Map 






FIG. 14: Top pane/; GRB positions as recorded by BATSE, 
with error bars indicated as light circles/ovals around the 
GRBs; pixelhzed at NSIDE=64 (note that the dynamic range 
is hmited to 2 even though a few pixels have 3 or 4 counts in 
them). Bottom panel: The BATSE exposure function, which 
is proportional to the amount of time spent monitoring a given 
direction and varies with declination, in Galactic coordinates. 



when gamma-ray count rates exceeded some minimum 
threshold relative to background in two or more of the 
eight detector modules, within some energy range [83] . 
Nominally, that energy range was 50-300 keV, but on the 
order of 30 percent of BATSE's observing time was spent 
with one of several trigger energy ranges different from 
this. In addition, while the minimum detection threshold 
in count rates relative to background was 5.5 sigma as a 
baseline, this value also changed many times over the 
course of the experiment, and was not always the same 
for different time intervals (BATSE tested count rates at 
64, 256, and 1024 ms intervals). 

A different trigger energy range essentially represents 
a distinct burst experiment, and a different detection 
threshold also changes the parameters of the experiment 
in an important way. However, we argue that the time 
variation in BATSE's ability to detect gamma-ray bursts 
is not sufficiently great to affect our results on the dipole, 
especially given the lack of statistics for the BATSE sam- 
ple of 2702 bursts. Kommers et al. [88] performed a 
search for gamma-ray bursts and other gamma-ray tran- 
sient phenomena with peak fluxes below (by a factor of 



~ 2) the flux necessary to count as a detection, and 
also with energies outside the nominal 50-300 keV. They 
found that the direction and intensity distributions of 91 
likely GRB candidates not included in the final BATSE 
catalog imply that biases associated with the trigger 
mechanism do not significantly affect the completeness 
of the catalog. In addition to this result, there is no rea- 
son to expect that changes in experimental parameters 
would have a particular effect on the dipole quantity we 
investigate here. We note especially that changing trig- 
ger criteria and energy ranges do not appreciably increase 
the chances that GRBs will be confused with Galactic 
sources (e.g., soft gamma repeaters (SGRs)), so contam- 
ination remains minimal; and any changes in trigger cri- 
teria apply uniformly over the entire sky, so there is no 
obvious reason why this would induce a dipole pattern. 
We therefore proceed to analyze the full catalog without 
accounting for these changes. 

However, one other experimental parameter is impor- 
tant for our dipole analysis: sky exposure in BATSE 
varies as a function of declination (BATSE spent differ- 
ent amounts of time looking at different declinations). 
We create a template out of the exposure function (see 
Fig. 14) and use this template as one of the systematic 



templates ti in the dipole formalism outlined in Sec. HI 



This corresponds to weighting pixels according to how 
much time the satellite spent observing a given area of 
the sky: see the approach in, for example, Scharf et al. 
[80] . The choices we make here, to ignore changes in 
trigger criteria but take the exposure function into ac- 
count, corresponds to the choices made in the paper by 
Tegmark et al. [81 calculating the angular power spec- 
trum of the BATSE 3B catalog, which found no evidence 
of deviations from isotropy on any angular scale. Our 
present tests can be regarded as updates (since we use the 
current catalog, which more than doubles the number of 
bursts in the 3B catalog) of Tegmark's, with focus on the 
dipole (direction as well as magnitude) using a real-space 
estimator that more naturally incorporates sky cuts and 
systematic templates. (Note that Tegmark et al. also im- 
pose a weighted averaging scheme in harmonic space to 
account for the very large position errors associated with 
GRBs, which are on the order of degrees, orders of mag- 
nitude larger than typical position errors associated with 
galaxies. This is unnecessary in our case given that the 
dipole probes scales much larger than the uncertainties 
in GRB positions.) 

The redshift selection function for GRBs is still only 
poorly understood, though better statistics are consis- 
tently being built up. GRBs come from even higher 
redshifts on average than NVSS sources (see Sec. VIII) 



(e.g., Xiao and Schaefer [89 ), however, so we can confi- 
dently say that regardless of the precise distribution, the 
local-structure dipole will be subdominant in compari- 
son with the kinematic dipole (given that it is already 
subdominant for the NVSS sources discussed in the next 
section), and so we consider only the kinematic dipole as 
a theoretical expectation below. 
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C. Systematic Checks and Dipole Amplitude 

The positions of GRBs detected by BATSE are shown 
in Fig.[l4] They are not very weh-locahzed; the positions 
typicahy have error bars on the order of degrees. The 
GRBs do not appear to cluster in any particular way by 
eye, but we apply our usual tests to see whether this 
holds up rigorously 

In considering what systematic templates to put in 
place, maps of Galactic foregrounds are unnecessary In 
particular, inclusion of the SFD dust template is unnec- 
essary since gamma rays are highly penetrating and not 
subject to appreciable dust extinction. We have explic- 
itly verified that the difference between the results includ- 
ing vs. not including SFD template is completely negligi- 
ble. In principle, there is the possibility of confusion with 
soft gamma repeaters (SGRs) or other sources of gamma 
rays (pulsars, terrestrial gamma-ray fiashes, black holes, 
etc.). This is highly unlikely given that GRBs are easy 
to distinguish from other gamma-ray sources based on 
spectral and time-domain data. However, any foreground 
objects that might contaminate a pure GRB sample are 
expected to vary with Galactic latitude, and since there 
is no reason not to do so, we run our usual test of progres- 
sively excising the Galactic plane. However, we expect 
no issues with astrophysical foregrounds given the rela- 
tively clean nature of GRBs as a source; in fact, no sky 
cut at all should be necessary. 

In all tests below, the systematic template we do use, 
as alluded to above, is the BATSE exposure function, 
which varies significantly with declination, mimicking 
(partially) a dipole (see Fig. 14). Some of the results 
presented below change both quantitatively and qualita- 
tively (e.g., statements we would make about the super- 
galactic plane are different) if this template is not taken 
into account. We also include quadrupole and octopole 
templates, as usual. This does not affect the results apart 
from widening the error bars. 

Given that cosmic variance should be much smaller on 
these scales than it was for 2MASS, we expect that re- 
sults at different values of minimum \b\ will be consistent 
with each other within the measurement error bars. Re- 
sults for varying Galactic 6cut are given in Fig. 15, from 
which it is clear that results at different 6cut are indeed 
consistent with each other. As expected, there is no de- 
tectable signal that varies with Galactic latitude. 

We now proceed to run the usual test cutting the su- 
pergalactic plane. It should be noted that for some cuts, 
a dipole is detected at marginal significance in this series 
of tests if the BATSE exposure function is not taken into 
account. Any kind of detection of the supergalactic plane 
in gamma-ray bursts would be very surprising given the 
complete lack of any association between GRBs and the 
local structure represented by the SGP (which goes out 
to something on the order of z ~ 0.02 or 0.03, depend- 
ing on estimates; see Lahav et al. [63]), but the result 
turns null when we account for the exposure function. 
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See Table VII for the results of cutting in supergalactic 



FIG. 15: Dipole amplitude in BATSE as a function of bcut- 
Note that all measurements of the dipole amplitude for differ- 
ent Galactic cuts are mutually consistent, and are consistent 
with zero as well. Measurement errors in the form of 68 and 
95 percent confidence intervals are shown; they tend to grow 
with the aggressiveness of the cut, as expected. 



latitude. 

We find another null result, which is good considering 
that detection of the supergalactic plane in GRB data 
would be a highly unusual find, and would almost cer- 
tainly indicate something problematic about our analy- 
sis. 



D. Dipole Direction and Conclusion 

The still-relatively-small number of gamma-ray bursts 
(2702) in the BATSE catalog places only the loosest of 
constraints on the direction of the dipole. All but a very 
tiny patch of sky (centered around (/,6) = (273°, 31°), 
and extending roughly 10 degrees in radius) is within 
the 3-sigma confidence interval for the direction of the 
dipole, and the CMB kinematic dipole direction (which 
is an appropriate direction of comparison in this case, 
since the kinematic dipole should dominate contributions 
to the dipole for GRBs) is only marginally outside the 2- 
sigma confidence interval (which given the looseness of 
the constraints is not a noteworthy result). 

In all, the BATSE data places useful constraints on the 
dipole amplitude, but not direction, in GRB data. Our 
results, while they are not yet strong enough to allow for 
the detection of the expected kinematic dipole, do place 
constraints on our ability to distinguish BATSE GRBs as 
a tracer of large-scale structure, and also constrain any 
intrinsic dipole in the large-scale structure that would 
manifest itself in the distribution of GRBs on our sky: at 
95 percent confidence, A < 0.117 for the intrinsic dipole. 
Because the constraints on the dipole are relatively weak 
relative to constraints using other surveys employed in 
this paper, and because the radial distribution of BATSE 
sources N{z) is poorly known, we have not attempted to 
produce a theoretical expectation for the BATSE local- 
structure dipole. 

In the course of this analysis we have, however, ef- 
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TABLE VII: Key patterns in cutting in supergalactic coordinates, for BATSE gamma-ray bursts. 
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fectively performed several sanity checks on the BATSE 
dataset, showing that the supergalactic plane is unde- 
tectable (as expected) using our analysis, that the Galaxy 
does not show up at all in the data (which is in line with 
previous studies of the GRB distribution), and that the 
BATSE exposure map must be taken into account in or- 
der for these tests not to turn up anomalous. We would 
of course still like to see constraints on the intrinsic dipole 
much better than those available from GRBs, those con- 
straints being on the order of 10~^. For this, we turn to 
the radio survey NVSS. 



VIII. DIPOLE IN NVSS 

For a long time, it was assumed that the distribution 
of radio sources was, like that of GRBs, indistinguish- 
able from isotropic and unclustered (e.g., Webster [90 ). 
In fact, even if the distribution of radio sources was not 
intrinsically isotropic, radio sources have a large range 
of intrinsic luminosities, and so structures would natu- 
rally wash out when sources were projected onto the sky 
and radial information was removed (Baleisis et al. [57 ). 
However, recent results, especially using the NRAO VLA 
Sky Survey (NVSS), have detected clustering in radio 
sources, and in particular, a dipole. NVSS is a radio sur- 
vey with nearly 1.8 million extragalactic sources (Condon 
et al. [9j). This survey presents an excellent opportunity 
to actually test for the presence of the kinematic dipole 
and possibly the intrinsic dipole in large-scale structure. 
NVSS has more potentially non-negligible systematics to 
control for than the other surveys we use, but it also 
has higher potential payoff because of its combination of 
depth and sky coverage. 



A. Previous Work 

Several attempts at calculating the dipole in radio 
sources have been made in recent years. 

Baleisis et al. [57^ present theoretical predictions and 
observational results for the dipole in the Green Bank 
1987 and Parkes-MIT-NRAO (PMN) catalogs. The com- 
bination of these catalogs gives ~ 40, 000 sources with 



flux > 50 mJy at 4.85 GHz. They find that the magni- 
tude of the dipole is an order of magnitude larger than 
expected from the contributions of large-scale structure 
(analogous to our local-structure dipole) and the kine- 
matic dipole. However, they are plagued by several sys- 
tematic errors. First, they find that the two catalogs 
they used have a mismatch in flux. While they correct 
for this, it is hard to do so with high precision and confi- 
dence. They also note that the radio sources in their cat- 
alogs are likely drawn from multiple populations, though 
this is true of any analysis that uses radio sources, and 
is not a crippling problem if the redshift distribution is 
sufficiently well-understood. 

Blake and Wall [56 attempt to measure the kinematic 
dipole alone in NVSS (see also Blake et al. for anal- 
ysis of the rest of the angular power spectrum in NVSS). 
They make efforts to remove the contribution of what 
they refer to as the "clustering dipole," the dipole that 
when flux- weighted gives a measure of the acceleration of 
the Local Group, and when unweighted matches up with 
our local-structure dipole. They claim that the cluster- 
ing dipole should die away by z < 0.03 (based on results 
from the Rowan- Robinson et al. [59^ analysis of the IRAS 
PSCz dipole, though the results of Erdogdu et al. [17] call 
this convergence into question) and contribute roughly 
A ^ 2 X 10~^ to the total amplitude of the total dipole, 
if it is not removed as they attempt to do. (Note that we 
have converted from their peak-to-trough "amplitude" S 
to our peak-to-zero amplitude A). 

Blake and Wall measure the remaining dipole - which 
would ideally be a kinematic dipole only, but which will 
in reality take contributions not only from local structure 
beyond what IRAS observed, but also from more distant 
large-scale structure, as we show below - by expanding 
the angular distribution of sources in spherical harmon- 
ics and measuring the harmonic coefficients airn for the 
dipole, quadrupole, and octopole, including all m values. 
Inclusion of higher harmonics is necessary because of the 
lack of full- sky coverage. They find that a dipole model 
is a good fit by and find good agreement with the 
direction of the CMB velocity dipole, which they cite as 
e = 97.2 ± 0.1°, = 168.0 ± 0.1°, and which converts to 
{l,b) = 264.3°, 48.1°). 

We take the Blake and Wall results as the most reli- 



28 



able previous result/^ and compare our results to theirs. 
However, we do have reason to expect that our estimate 
of the dipole will be more reliable than theirs; in par- 
ticular, we do a great deal more to take systematic ef- 
fects into account than they do in their analysis, using 
our real-space estimator. We also do not attempt to re- 
move local sources as they do, since our objective is to 
compare our observational results to the full dipole sig- 
nal expected from theory, which includes both a local- 
structure and a kinematic contribution. Since we are not 
flux-weighting, local sources do not contribute preferen- 
tially to the dipole, and so we can afford to leave them 
in the analysis. 



B. Theoretical Predictions 

For NVSS, we must be careful to include in our theo- 
retical predictions the contributions of not only the local- 
structure dipole, but also the kinematic dipole. It is no 
longer the case, as it was for 2MASS and 2MRS, that the 
kinematic dipole is swamped by two orders of magnitude 
by the local-structure dipole. Rather, the two are on the 
same order of magnitude, as also recognized by Baleisis 
et al. [57] and Blake and Wah [56 . 

The local-structure contribution is calculated in the 
same way as always. This part of the prediction will 
vary somewhat depending on what redshift distribution 
n{z) we use. Dunlop and Peacock [94 derived n{z) for 
several different flux cutoffs, though the results for the 
dipole amplitude and the redshift distribution itself are 
somewhat robust to changes in the flux cutoff since radio 
galaxies display such a wide range of intrinsic luminosi- 
ties (Baleisis et al. [57 ; Blake et al. [92] note specifically 
that for the NVSS frequency of 1.4 GHz, the clustering 
of radio galaxies is not strongly dependent on flux for 
fluxes between 3 mJy and 50 mJy). The redshift distri- 
bution developed by Ho et al. [31] as a best model for 
NVSS avoids several drawbacks of the Dunlop and Pea- 
cock distribution, especially the assumption that bias is 
redshift-independent and the heavy reliance on the func- 
tional form of the luminosity function rather than data in 
constraining the redshift distribution. However, without 
repeating Ho et al.'s rather detailed analysis, we are left 
unable to calculate the redshift distribution for flux cuts 
different than the 2.5 mJy cut that they use. 



Fortunately, we can here make use of the fact that the 
predictions are not strongly tied to the flux cutoff. We 
therefore follow Ho et al. in modeling the NVSS redshift 
distribution as follows: 



W{z) 



z I 



(33) 



where W{^z) = 6(z)n(z), where h{^z) is the bias as a func- 
tion of redshift and n{z) (H(z) in Ho et al.) is the prob- 
ability distribution for the galaxy redshift. Ho et al. give 
6eff = 1.98, = 0.79, and a = 1.18 as best-fit parame- 
ters. Using all of this, we find that for this distribution, 
the contribution of the local-structure dipole to the total 
dipole in NVSS is A = 9.8 x 10""^ - 0.0010. 

Meanwhile, the theoretically expected kinematic 
dipole may be calculated as shown in Appendix [Aj The 
predicted amplitude is 



A = 2^ = 2[1 + 1.25x(l - p)]P 



(34) 



where the first term in brackets essentially represents 
the contribution of relativistic aberration and the second 
term represents the contribution of the Doppler effect. 
From the CMB, (3 = v/c = 1.23 x 10"^. Meanwhile, 
X and p (exponents in the power laws for the intrinsic 
number counts of galaxies as a function of limiting mag- 
nitude, and for the intrinsic flux density of a galaxy as 
a function of frequency, respectively; see Sec. HC) are 



not known precisely, but can be estimated, as we do just 
below. 

This expression is equivalent to that used by Blake and 
Wall [56] (converting from their S to our A = S/2): 



A= [2 + x(l + a)]/3 



(35) 



with the substitution in Eq. ( |A9[ ) x 2x/b and p —a. 
(The latter substitution is a straightforward matter of no- 
tation; the first has to do with switching from magnitudes 
to fluxes.^2) We follow Blake and Wah (for NVSS) and 
also Baleisis et al. (who were not working with NVSS, but 
who did work with radio catalogs) and take x ^ 1 and 
a ~ 0.75, which yields a kinematic dipole amplitude of 
Akin ^ 0.0046. The full theoretical value includes this 
kinematic dipole and the contribution from the local- 
structure dipole (Aiocstr = 0.0010). We can therefore 
write the full theoretical prediction as 



A = 0.0046 ± 0.0029 ± 0.0006, 



(36) 



Singal |93| has performed the most recent analysis in this vein. 
His results are suspicious, as he finds truly exorbitant speeds for 
the Local Group (on the order of 1700 km/s). This becomes 
more understandable given that the way in which he accounts 
for the sky cut, particularly the hole in NVSS at declination 
< 40°, is suspect (he simply cuts out dec > 40° as well in or- 
der to counterbalance the hole at < 40°). Also, his method of 
detecting the dipole does not account for coupling between the 
dipole and other multipoles on the cut sky, and also neglects any 
contribution from the local-structure dipole to the results. 



^ In our original notation, which follows Itoh et al. the number 
of galaxies detected by the survey is proportional to 10^"^ii™, 
where miim is the limiting magnitude of the survey. In Blake and 
Wall 56 , the number of galaxies is proportional to , where 
S is flux, and we have added a prime to distinguish variables 
named x. Equating with lO^^ii"^ (up to a proportionality 

constant) and using the fact that S oc 10~^^ii"^/^, we have that 
2m\iuiX /5 = xmiim, or x = 2x' /5. Blake and Wall give x ^ 1. 
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NVSS, All Sources 




NVSS, Sources With Flux Greater Than 15 mJy 




FIG. 16: Top panel: All sources in NVSS, in Galactic coordi- 
nates. Note the "hole" in the data for declinations less than 
40 degrees, and the declination-dependent striping (visible in 
this coordinate system as a series of "wavy" stripes going out- 
ward from the pattern set by the declination-dependent hole 
in the data). Bottom panel: Sources with flux greater than 
15 mJy. The spurious power goes largely away with this flux 
cut, and can no longer be seen by eye. (Dynamic ranges are 
restricted in both these maps so as to better show structure.) 



where the first uncertainty corresponds to the kinematic 
dipole and the second to the local-structure dipole. If we 
find a result that is outside these cosmic-variance errors 
from the central value, and this is not a systematic effect, 
we might invoke the presence of an intrinsic dipole as an 
explanation. 



C. Present Work 

We turn first to examining the systematics that need 
to be accounted for in the NVSS data. These systematics 



are illustrated in Fig. 16 

First, the survey did not observe below a declination 
of 40 degrees. For our purposes, the pixellization around 
this "hole" in declination is especially important to pay 
attention to, as we find that a sky cut for pixels with 
centers at dec < 40° gives significantly different results 
for the dipole amplitude and direction than a cut at dec 
< 37° at low-resolution pixellization (NSIDE=16). We 
work at much higher resolution (NSIDE=128), where the 
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FIG. 17: Measured NVSS dipole amplitude as a function of 
flux cut, with error bars. We show both our results and those 
from Blake and Wall [56 . For the dotted red curve, we include 
the Haslam et al. map of 408 MHz Galactic emission as a sys- 
tematics template, and remove \b\ < 15° and dec greater than 
78 degrees/less than -10 degrees; for the solid green curve, we 
include sources with dec greater than -37 degrees (thus ef- 
fectively increasing /sky from 0.42 to 0.60), but also include 
quadrupole, octopole, and dec-dependent-striping templates. 
We suspect that the near-monotonic increase of the ampli- 
tude with the flux cut is due to spurious power in the NVSS 
map, and possibly the presence of local structure in the sur- 
vey as well (since local structures preferentially have higher 
fluxes). The apparent agreement between theoretical predic- 
tions and the Blake and Wall results is partially misleading, 
as discussed in the text. 



effect is not as strong, but we still cut all pixels with dec 
< 37° just to be safe. We choose that particular number 
following Smith and Huterer [95 , and also because this 
appears to be a good conservative choice if we want to 
completely avoid problems associated with the hole in 
declination. 

Second, faint sources very close to bright sources can- 
not be reliably detected. We therefore mask out a 0.6° 
radius around all sources brighter than 2500 mJy, fol- 
lowing Ho et al. [31]. Blake and Wall do not perform 
this same masking, although they do remove known lo- 
cal sources. 

Third, we use as a systematic template the map of 
Galactic synchrotron radiation of Haslam et al. [%] and 
Haslam et al. [97 . This is a 408 MHz radio continuum 
map of the entire sky that combines data from four dif- 
ferent surveys and is dominated by synchrotron emis- 
sion. Information about foreground emission is impor- 
tant when dealing with a radio survey such as NVSS, and 
we can apply this map as a systematic template since 
NVSS could plausibly pick up non-extragalactic signal 
from this emission. 

Fourth, there is declination-dependent "striping" in 
the maps, which appears even by eye if no flux cut is im- 
posed. This problem stems from the fact that the NVSS 
observations were made using two different configurations 
of the VLA, the D configuration for observations between 
declinations of —10° and +78°, and the DnC configura- 
tion for declinations between —40° and —10°, and above 
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+78°. The striping is readily apparent by eye for the fuh 
catalog, but is invisible by eye and largely absent even 
in more rigorous tests for fluxes above ~ 15 mJy (see 
Blake and Wall [56]). We therefore begin by examining 
the stability of our dipole results as a function of flux cut. 

We find that neither the direction nor the amplitude 
of the dipole is stable for different flux cuts. This is 
as expected for flux thresholds less than 15 mJy, since 
the striping artifact gradually dies away as the lower 
flux threshold is increased from zero up to 15 mJy. 
But for flux thresholds above 15 mJy, the fact that the 
dipole remains unstable is a problem. Increases in the 
dipole amplitude might be due to the influence of the 
local- structure dipole, which could come into play more 
strongly for these brighter (and therefore at least some- 
what more local) sources. Fluctuations in the direction 
(the best-fit / ranges from 219.1 to 234.3, and the best-fit 
b ranges from 11.9 down to -0.2) could also be the result 
of the local-structure dipole, though this seems unlikely 
since the b coordinate in particular moves away from the 
direction of the local-structure dipole that we have seen 
in previous tests using 2MASS as we go to brighter and 
brighter sources (higher and higher flux thresholds) in 
NVSS. 

We implement two distinct and complementary strate- 
gies for dealing with these issues. 

(1) The first strategy is to make a more aggressive cut 
in declination, which makes use of only one subset of 
the NVSS maps, corresponding to the D configuration of 
the VLA: that is, we remove all portions of the sky with 
declination less than -10° or greater than 78°. These 
results show a great deal more stability in the direction of 
the dipole for flux cuts above 5 mJy, and certainly above 
15 mJy. The primary results can be seen in Fig. We 
have also examined the stability of results as a function 
of cut in Galactic b for the fixed case of a flux cut at 25 
mJy. This provides our usual test for contamination that 
varies as a function of Galactic latitude, and helps justify 
our choice of \b\ < 15° for our cut. We find that there 
is some fluctuation in the dipole signal as a function of 
6cut5 likely an indication of remaining spurious power in 
the map rather than genuine Galactic latitude-dependent 
contamination. 

We regard the result with a flux cut at 15 mJy and 
\b\ < 15° as paradigmatic, since it provides the best com- 
promise between a large number of sources (getting more 
sources requires getting into flux ranges where the re- 
sults are less trustworthy) and having stable results for 
the dipole. Using this result as the fiducial, we perform 
our usual test of making cuts in supergalactic latitude, 
and find that the presence of the supergalactic plane is 
not visible in the results: the results are extremely stable 
as a function of cut in SGB, with the ratio /sky//sources 
changing by less than 1 percent in all but one of the SGB 
cuts we study. We also compute our likelihood distribu- 
tion P{Cl^\Cf'^) using the combination of a 15 mJy flux 
cut and \b\ < 15° cut. We calculate the direction of the 
dipole for this same case, compare with the CMB kine- 
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FIG. 18: Top panel: Likelihood of dipole direction in NVSS, 
marginalized over amplitude. Bottom panel: The CMB dipole 
direction is over 3 sigma away from the best-fit NVSS dipole 
direction. This is not a problem, however, since the NVSS 
dipole amplitude is A ^ 10~^, and we expect the LSS dipole 
direction to match the CMB dipole direction only when ampli- 
tudes A ^ 10~^ are probed. Both of these plots correspond to 
the case in which sources have flux greater than 15 mJy, only 
declinations between —10° and 78° are kept, and \b\ < 15° is 
cut. 



matic dipole, and show the results in Fig. 18 The final 
plot of this section is Fig. [19) which depicts the different 
cuts used in various portions of the analysis for NVSS. 

(2) The second strategy for attempting to control the 
various systematics in NVSS is to include systematics 
templates that consist of degree-wide stripes in declina- 
tion: that is, rather than making highly selective cuts in 
declination as before, we instead keep the entire NVSS 
map above —37° declination and then include system- 
atics templates that have value 1 within the declination 
range (l9, 6> + 1], 6> = -37°, -36°, 88°, 89°, and are zero 
elsewhere; this is a total of 127 additional templates. 
This allows us to keep a much greater percentage of the 
sky in the analysis, and correspondingly more sources. 
This also allows us to marginalize over the quadrupole 
and octopole modes by incorporating them as systemat- 
ics templates too (we avoided doing this with the more 
aggressive cut in declination since using so many tem- 
plates on such a small fraction of the sky inflates the 
dipole amplitudes and error bars). The results of this 
test are also shown in Fig. [TTl 

In this case, the single most reliable value for the dipole 
amplitude is A = 0.0280, the value obtained from analyz- 
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TABLE VIII: Summary of most reliable single results from each survey. From left to right in the table appear the name 
of the survey, the redshift range probed by the survey, the fraction of the sky covered, the number of sources available in 
the most reliable subset of the dataset, the observed dipole amplitude with error bar, the theoretical dipole amplitude (with 
cosmic- variance error bar if applicable), the direction of the best-fit observed dipole in Galactic coordinates (/, 6), and the most 
important systematic effect (in some cases, out of several candidates) that must be taken into account in attempting to detect 
a dipole in the dataset. 



Survey 


Redshift 


/sky 


N 


^obs 




{I, b)ohs 


2MRS 


< z < 0.1 


0.86 


41,834 


0.120 ±0.009 


0.311 ±0.122 


(213.8°, 35.2°) 


2MASS 


0<^<0.2 


0.65 


386,008 


0.104 ±0.004 


0.084 ± 0.033 


(268.4°, 0.0°) 


BATSE 


^>2 


1.00 


2702 


< 0.051 (68% CL) 


[unc. prediction] 


[weak constraints] 


NVSS 


z ^ 1 


0.42 


211,487 


0.027 ±0.005 


0.0046 ± 0.0035 


(214.5°, 15.6°) 



ing the sample with fiux greater than 2.5 mJy. While the 
dipole amplitude changes as a function of flux cut, the 
direction is stable when we keep sources above —37° dec- 
lination and use templates to correct for the declination- 
dependent striping, and we no longer have any clear rea- 
son to doubt the results drawn from the most permissive 
flux cuts. The major caveats in this case are that (a) the 
use of "stripe" templates may effectively eliminate actual 
signal as well as spurious power, and (b) the direction we 
obtain for the dipole in case (2) does not agree with that 
from case (1). 

It is, however, at least somewhat encouraging that 
the "best-case" observed amplitude of the NVSS dipole 
is similar in both cases: A = 0.0272 in case (1) and 
A = 0.0280 in case (2). The agreement between these 
values may or may not be coincidental, but either way, 
the value is not consistent with the theoretical prediction 
of A = 0.0046, at well over 99 percent confidence. The 
reasonable implication is that the dipole we are measur- 
ing in the NVSS map is partially spurious. There are 
several possible sources of this spurious signal: the two 
most likely are that the declination-dependent striping 
may not be wholly taken care of in the maps we test, 
even with our aggressive flux cuts in case (1) and our 
templates in case (2) (observational error); or that local 
sources contribute more strongly to the local-structure 
dipole than our theoretical modeling allowed for, thus 
lifting the local-structure dipole's contribution to the or- 
der of 10~^ rather than 10~^ (theoretical error). 

We strongly suspect that observational error plays the 
much greater role. Theoretical error seems like a reason- 
able explanation in case (1), given that the theoretical 
prediction for the dipole amplitude was computed us- 
ing a flux cut of 2.5 mJy, while the observational results 
were computed in case (1) using a flux cut of 15.0 mJy. 
But it does not explain the discrepancy in case (2), and 
in any case the theoretical predictions are ideally fairly 
robust to flux cuts. (That said, for the most restric- 
tive fiux cuts, there is an order-of-magnitude difference 



in fiux between theory and observation, and this is likely 
to motivate adjusting the theoretical predictions at least 
somewhat, adding back in the local-structure dipole that 
we expected to be subdominant.) 



As a final note, we also include in Fig. [T7| the results 
that Blake and Wall [56] found when doing a related 
dipole analysis on NVSS. There are a handful of reasons 
why the lack of agreement between our results and theirs 
is not unexpected. First, they remove local structure in 
an effort to search for the kinematic dipole only; since 
we do not do this, our results need not recover theirs, 
and in fact are expected to give a higher signal. Also, 
our cut in declination is different than theirs, leading to 
further potential discrepancies. The apparent agreement 
between theoretical predictions and the Blake and Wall 
results is also partially misleading, in that the theoreti- 
cal prediction includes contributions from both the kine- 
matic dipole and the local-structure dipole at redshifts 
not excised by Blake and Wall {z > 0.03), but they were 
attempting to measure only the kinematic dipole. 

All this said, it is clear that radio surveys of this sort 
are an excellent setting for the tests we perform, and we 
look forward to maps from, e.g., LOFAR and SKA to 
perform similar tests. 



IX. CONCLUSION 

In this work, we have focused on what might be called 
the most straightforward tests of statistical isotropy in 
large-scale structure - looking for dipole signals in exist- 
ing surveys over a wide range of wavelengths. It turns 
out that, despite the relative straightforwardness of the 
tests themselves (Sec. Ill), the results must be carefully 
interpreted, as dipole signals take contributions from sev- 
eral different sources. Some of these sources, such as 
local-structure and kinematic dipoles, are theoretically 
quite well-understood, while others, such as the intrinsic 
dipole, may be less so; see Sec. Ill) 
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\b\ < 15°; \b\ <30° 




same as immediately above but with 
|SGB| > 41.15° also cut 



FIG. 19: Various cuts employed in the analysis of the NVSS 
map. The background map is the NVSS map with all the 
sources that have flux greater than 15 mJy. 



Observational results in infrared (2MRS/2MASS; 
Sec.|lV|and|vl), gamm a rays (CGRO/BATSE; Sec.[VIl|), 
and radio (NVSS; Sec. VIII) turn up no seriously unex- 



pected results, in either dipole amplitude or dipole di- 
rection. As long as we are careful to take all sources of 
dipole signal into account in our theoretical modeling, 
the observational results are in line with theoretical pre- 
dictions. 

Rigorous tests of this sort, while they ultimately turn 
up no unexpected results, are valuable tests of current 
cosmological models, as they add new wavebands in 
which rigorous tests of statistical isotropy have been con- 
ducted, and ensure that statistical isotropy is probed at 
different epochs, using different surveys with different 
systematics. Combined with similar tests using maps of 



the cosmic microwave background, these measurements 
impose interesting constraints on models of, and physi- 
cal processes during, cosmological inflation. We provide 
a summary of the most basic results, which are elabo- 
rated upon heavily in the body of the paper, in Table 



One feature of these results is particularly worth high- 
lighting: namely, that they place constraints on the am- 
plitude of any intrinsic dipole present in large-scale struc- 
ture. This is especially true of the BATSE and NVSS re- 
sults, since they are not expected to be dominated by the 
local-structure dipole. BATSE places an upper bound 
on the intrinsic dipole amplitude at 1 x 10~^ at 95 per- 
cent confidence, while our most optimistic NVSS results 
places an upper bound at 4 x 10~^ at 95 percent confi- 
dence. As discussed in Sec. VIII C[ however, we cannot 
place as much confidence in the NVSS result as we do 
in results from the other three surveys analyzed in this 
paper. Detection of a dipole signal in any survey requires 
that a great deal of attention be devoted to controlling 
for systematic errors and spurious power in the survey, 
and while we have considerable confidence that we have 
done this success fully for 2MRS, 2MASS, and BATSE, 
we outline in Sec. |VIIl"C] why acquiring and interpreting 
NVSS results presents the greatest challenge of all. 



Other surveys may provide interesting candidates for 
these same kinds of tests in the future. Sloan luminous 
red galaxies (LRGs) are a sufficiently clean dataset that 
these tests may be applicable and workable there (see, 
e.g.. Abate and Feldman [98 for a related test). Also, 
survey results that have yet to be released may be use- 
ful. The Wide-field Infrared Survey Explorer (WISE), 
which produced a preliminary data release in April 2011 
and has full-survey results forthcoming now, covers more 
than 99 percent of the sky (Wright et al. [99 ). Tests per- 
formed on WISE would be similar to tests performed on 
2MASS, but would update 2MASS results with a more 
recent and much deeper survey. As the X-ray background 
becomes better-understood, this may also serve as an in- 
creasingly valuable test of statistical isotropy and set- 
ting in which to attempt to detect dipole signals. The 
Dark Energy Survey (DES) will be useful in probing the 
distribution of galaxies to high redshift, and will have 
sufficient sky coverage to make the tests presented here 
useful. In microwaves, dipole signals might be detectable 
in maps of the gravitational lensing of the CMB, which 
provide a very good tracer of mass. Finally, new radio 
surveys such as LOFAR and SKA will probe orders of 
magnitude more sources down to far lower flux thresh- 
olds than NVSS (see, e.g., Crawford |lQQj ). and would 
provide very valuable updates to NVSS results on the 
dipole amplitude and direction. The kinematic dipole, 
both direction and amplitude, should be unambiguously 
detected in these surveys. 



33 



Acknowledgments 

For help with the set of research projects that have 
gone into this paper, thanks to Jim Condon, Sudeep 
Das, Jerry Fishman, Chris Hirata, Michael Hudson, Tom 
Jarrett, Lucas Macri, Dominik Schwarz, and Weikang 
Zheng, all of whom were very helpful at different points 
along the way. For useful comments in the aftermath of 
initial arXiv posting, we thank Maciej Bilicki and Jim 
Zibin. CG would additionally like to thank the members 
of his dissertation committee, all of whom provided valu- 
able input: Philip Hughes, Jeff McMahon, Aaron Pierce, 
and especially Tim McKay. DH would additionally like 
to thank Glenn Starkman for sparking his interest in test- 
ing statistical isotropy with large-scale structure. 

This publication makes use of data products from the 
Two Micron All Sky Survey, which is a joint project of the 
University of Massachusetts and the Infrared Processing 
and Analysis Center/California Institute of Technology, 
funded by the National Aeronautics and Space Admin- 
istration and the National Science Foundation. We also 
acknowledge the use of the publicly available HEALPix 
[54] package. 

We have been supported by DOE grant under con- 
tract DE-FG02-95ER40899, NSF under contract AST- 
0807564, and NASA under contract NNX09AC89G. DH 
thanks the Aspen Center for Physics, which is supported 
by NSF Grant No. 1066293, for hospitality. 



Appendix A: Relativistic Aberration and Doppler 
Effect 



We first address aberration, following the formalism of 
Buries and Rappaport |101j (see also, e.g., Calvao et al. 
[102 ), who derive equations for aberration with the ulti- 
mate goal of showing that aberration of the CMB temper- 
ature might be detectable statistically by Planck, looking 
at shifts of CMB peaks. While this is not our goal, the 
formalism still holds. 

We define a spherical-coordinate system with the z- 
axis in the direction of motion. If we take the "unprimed" 
frame to be the CMB frame, and the "primed" frame to 
be the frame of the Solar System barycenter, then the 
azimuthal angle (\) is unchanged between frames: (j) = (j)' . 
However, the polar angle is affected as follows: 



smc 



smf 



7(1 -/3 cos i90 



(Al) 



where j3 is the relative velocity of the Solar System with 
respect to the CMB, 7 = (1-^^)"^/^ as usual, and 6> = 
corresponds to the direction of forward motion. With the 
assumption that P is small, which is a good assumption 
given that CMB observations show it to be on the order 
of 10~^, expansion in a Taylor series gives 



Finally, we take the arcsin of both sides and expand the 
arcsin function assuming small /3 to obtain 



l9 = 6>' + /3 sin l9^ 



(A3) 



We are ultimately interested in calculating how areas 
(and volumes) on the celestial sphere are stretched or 
compressed, and hence want the quantity d£t/d£t' . With 
that in mind, we compute 



sm 



1 + /3COS 0' 
1 + /3COS 0' 



(A4) 
(A5) 



and find that dQ/dQ^ = (1 + PcosO'Y . Hence areas and 
volumes on the sky, proportional to sm{6)d6d(j), change 
as (1 + /3 cos i90^ ~ 1 + 2/3 cos i9'. 

Itoh et al. [11] provide a more complete derivation of 
this, including both the Doppler effect (which changes 
frequencies and hence measured magnitudes since we 
never measure bolometric magnitudes) and relativistic 
aberration, and derive the following expression for the 
observed angular number density of galaxies n{6) given 
the limiting magnitude miim^ 



m < miim) = n{m < mum) 



1 + 2/3 cos a 



where 



P= [l + 1.25x(l-p)]/3. 



(A6) 



(A7) 



Here the intrinsic flux density of a galaxy is assumed to 
be a power law 6'rest(^) ^ and the intrinsic number 
counts of galaxies n is 



n{m < m\im) oc 10^ 



(A8) 



where x is a numerical coefficient of order unity. The 
angle a is the angle between the angular direction 6 and 
the angular direction of the Earth's peculiar velocity v 
on the celestial sphere, the same as 0\ but with more 
convenient notation. The factor of 2 in 2^5 above comes 
from the same source as the squa re in (1 + /3cos6>')^ ear- 



lier. The correction for P in Eq. (A7) is the contribution 



of the Doppler effect to the overall kinematic dipole in 
observations of LSS. 

The final formula for the combined effects of rela- 
tivistic aberration and the Doppler effect is dQ/dQ' = 
1 + 2/3 cos 6\ so that the predicted amplitude of the kine- 
matic dipole is 



A = 2^ = 2[l + 1.25x(l-p)]/3. 



(A9) 



Appendix B: Quadrupole and Octopole Templates 



sin (9 = sin i9'(l + /3 cos i90. 



(A2) 



In Eq. ( 11 ) and below, we show how the dipole formal- 
ism from Hirata [7 can be used to separate out genuine 
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dipole signal in a map of large-scale structure from spuri- 
ous signal due to systematic effects or coupling with the 
monopole. In this Appendix, we detail how we can guard 
against the possibility that some of the "genuine" dipole 
signal actually comes from higher multipoles. 

The most straightforward way of doing this is simply 
to include the five 1 = 2 and seven ^ = 3 modes as sys- 
tematics templates in the analysis, corresponding to the 



^^/cit^(n) term in Eq. (11). More specifically, the tem- 



plates ti{u) are assigned^to be l2m(n) for — 2 < m < 2 
and l3^(n) for — 3 < m < 3. We expect that quadrupole 
and octopole modes should not contribute to the dipole 
signal at all in the limit of full- sky coverage. However, 
as the sky is cut, multipoles become coupled, and in the 
limit of very small sky coverage, there is high degener- 
acy between the quadrupole and octopole modes and the 
dipole mode, which is the only one of interest. 

To quantify this precisely, we run a simple test, the 
results of which are shown in Fig. 20 We take an 



artificially-generated map that contains nothing but a 
pure dipole in a certain direction (in this case, / = 
61.4°, 6 = 33.4°). The dipole amplitude is A = 0.1, 
so we expect that in the full-sky limit, our dipole es- 



timator should recover the result of A = 0.1, and the 
error bars should be the same regardless of whether we 
include quadrupole and octopole templates - there is, af- 
ter all, no coupling between the dipole and ^ = 2/^ = 3 
modes. However, as we make more and more aggressive 
cuts (accomplished here by making cuts that are symmet- 
ric in Galactic latitude, as we frequently do throughout 
our analyses of real datasets), we find that including the 
quadrupole and octopole templates becomes more and 
more important. While all results are consistent with 
the correct amplitude of A = 0.1, the error bars are much 
larger when i = 2 and i = 3 templates are included when 
roughly half the sky or more is cut. 

Given that the results depend only very weakly on in- 
clusion of quadrupole and octopole templates for small 
sky cuts, we do not always incorporate quadrupole and 
octopole templates into our analysis in the limit of nearly 
full- sky coverage. Meanwhile, especially in the case of 
NVSS where our analysis deals almost exclusively with 
less than half the sky, we sometimes explicitly compare 
the case where quadrupole and octopole templates are 
not included to the case where they are, since inclusion 
of the templates substantially weakens our results. 
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FIG. 20: Sensitivity of dipole amplitude measurements on 
marginalization over the quadrupole and octopole compo- 
nents, as a function of the fraction of the sky covered. We 
show the dependence of the detected dipole amplitude on /sky 
for the ideal case where the map being analyzed consists of a 
pure dipole. The dipole has amplitude A = 0.10, indicated by 
the dotted magenta line. We find that while the correct value 
of the amplitude is, in all cases, recovered within the appro- 
priate error bars, the inclusion of quadrupole and octopole 
templates does affect the size of the error bars. When /sky is 
more than about 0.7, the quadrupole and octopole templates 
make little difference, but below that value of /sky, error bars 
start increasing noticeably as more of the sky is cut away. 
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